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METHOD AND APPARATUS FOR SELECTION AND EVALUATION OF SUBSTANCES IN TREATMENT OF 

BONE DISORDERS 

5 

The present invention relates in general to different aspects of a more precise determina- 
tion of which types of defects or conditions may take place or exist in bones and how to 
determine which substance or substances to administer in order to counteract such de- 
fects or conditions. A large group of such defects or conditions is seen as different types 
10 of osteoporosis. 

Pharmacological agents used to manage osteoporosis can be divided according to their 
mechanism of action. Today two main groups exist, which are the antiresorptive agents 
and the bone formatting agents. Beside these two main groups there is a third group with 

15 agents having more heterogeneous effects such as the anabolic steroids, vitamin D and 
its analogues. Most of the agents that are commercial available today belongs to the 
antiresorptive group, which are believed to act by decreasing the rate of bone resorption 
and thereby slowing the rate of bone loss. Some of the agents have in large clinical stu- 
dies shown to reduce fracture risk. The optimal timing and duration of treatment for most 

20 of the agents remain to be defined and will probably differ between them. 

Below the different pharmacological agents, who may have a potential in relation to the 
management of osteoporosis, are grouped according to their mechanism of action. Only 
very few of the agents are approved for treatment of osteoporosis and most of them are 
25 still under development. Examples from each group are given together with the names of 
some of the producers. 

Antiresorptive Agents 

30 The antiresorptive agents are expected to act through a direct or indirect inhibition of the 
osteoclasts. The group may be divided into the following four classes: 
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Bisphosphonates 

Etidronate (Proctor & Gamble, Gentili) 

Clodronate (Astra, Boehringer Mannheim/Roche, Rhone-Poulenc Rorer) 
5 Pamidronate (Novartis, Gador) 

Alendronate (MSD, Gentili) 

Risedronate (Proctor & Gamble, Takeda) 

Ibandronate (Boehringer Mannheim/Roche) 

Zoledronate (Novartis) 
10 Tiludronate (Sanofi) 

Incadronate (Yamanouchi) 

Neridronate (Gentili) 

Olpadronate (Gador) 

EB-1053 (Leo) 
15 YH 529 (Yamanouchi - Hoechst Marion Roussel) 

Hormone Replacement Therapy 

Estrogens (Novartis, Novo Nordisk, Organon and others) 
20 Estrogens with progestogens (Novartis, Novo Nordisk, Organon and others) 
Tibolone (Organon) 

Trimegestone (Hoechst Marion Roussel) 
Selective Estrogen Receptor Modulator (SERM) 

25 

Raloxifene (Eli Lilly) 
Droloxifene (Pfizer) 
Idoxifene (SmithKline Beecham) 
Levormeloxifene (Novo Nordisk) 
30 Tamoxifen (Zeneca, Rhone-Poulenc Rorer and others) 

Calcitonin 

Salmon calcitonin (Novartis, Rhone-Poulenc Rorer, Leo and others) 
35 Human calcitonin (Suntory) 
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Calcitonin gene-related peptide (Warner-Lambert) 

The agents that so far have been approved for treatment of osteoporosis are nearly all 
belonging to the group of antiresorptive agents. 

5 

Agents with Heterogeneous Effect 

The agents in this group have so far been less well characterised with respect to their 
mechanism of action on the bone. 

10 

Stanozolol 

Oxandrolone 

Nandrolone 

Dihydrotachysterol and other vitamin D derivatives 
15 Alfacalcidiol 
Calciferol 
Calcitriol (Roche) 
Cholecalciferol 

ST-630 (Discovery Laboratories) 
20 Calcium (Novartis and others) 

Calcium with vitamin D (Nycomed and others) 
Ipriflavone and other isoflavones 
Strontium (Servier) 
Osteoprotegerin (OPG) (Amgen) 

25 

The Bone Formatting Agents 

The bone formatting agents are expected to act through a direct or indirect stimulation of 
the osteoblasts through a mitogenic mechanism. So far only the fluorides have exten- 
30 sively been used in the treatment of osteoporosis. Examples of bone formatting agents 
are given below: 

Sodium fluoride (Mission Pharmacal) 
Monofluorophosphate 
35 Parathyroid hormone^ (Eli Lilly) 
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Parathyroid hormonet^Allelix Biopharmaceuticals) 
Growth hormone (Novo Nordisk, Eli Lilly, Pharmacia & Upjohn and others) 
Growth hormone releasing compounds (Novo Nordisk and others) 
IGF-1 (Pharmacia & Upjohn, Novartis) 
5 IGF-BP3 (Celtrix Pharmaceuticals) 

Osteogenic potein-1 (Creative BioMolecules) 

Bone morphogenetic protein-2 (Genetics Institute, Yamanouchi) 

Because osteoporosis is associated with low bone mass, traditional bone mineral density 
10 (BMD) measured by Dual-energy X-ray absorptiometry (DXA) are used to assess the effi- 
cacy of treatment. However, there is growing evidence that BMD is just one of many fac- 
tors that contribute to the bone strength. Clinical studies have demonstrated that there is 
no simple relationship between changes in BMD induced by treatment and changes in 
fracture incidence. In studies with sodium fluoride, BMD at the lumbar spine increased by 
15 35% but there was no decrease in fracture incidence. Bisphosphonates have been shown 
to increase bone mass by about 10 %, which was associated with a reduction in fracture 
incidence of 50%. However, on the basis of a simple mass-strength correlation an in- 
crease of 10% in BMD after bisphosphonate treatment should only result in an increase in 
strength of 15% and, presumably, a corresponding decrease in fracture incidence. Also in 
20 relation hormone replacement therapy (HRT) and calcitonin significant reductions in frac- 
ture incidence have been demonstrated despite relatively small changes in BMD. The re- 
sults from several clinical studies have shown that fracture incidence can be reduced sig- 
nificantly with only a small increase in BMD, and that this discrepancy may be explained 
by an effect of treatment on the bone micro-architecture. 

25 

Bone micro-architecture can be evaluated though histomorphometry based on biopsies 
drilled out of the ilium bone. This method is very costly and time consuming and also in- 
convenient to the patients and therefore only applicable in small experimental studies. It 
has been suggested that quantitative computed topography (QCT) could be used to as- 

30 sess changes in the bone quality, but due to the high precision error with this technique it 
is presently not usable for follow-up measurements. Further the QCT devices are very ex- 
pensive and also the higher x-ray dose limits the use for routine measurements. Besides, 
looking at the fracture incidence, the only non-invasive technique available today in rela- 
tion to assessing the effect of an anti-osteoporotic treatment is by measuring BMD. A 

35 newly developed method termed Digital X-ray Radiogrammetry (DXR) provides more 
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relevant and detailed information about treatment effect than can be obtained with con- 
ventional DXA technology. Data from a recent pilot study seems to indicate that the 
qualitative bone parameters generated by DXR is more sensitive than DXA to changes in 
the bones following anti resorptive treatment. Further, the DXR technology is also ex- 
5 pected to be able to distinguish between the individual anti resorptive agents with respect 
to their pharmacological effects on the bones. 

Only a few investigations have been performed in order to determine the actual effect of 
such substances on specific parameters. Normally, these investigations entail firstly pro- 

10 viding a certain condition in a large number of test animals and subsequently administer- 
ing the substance. The actual determination of the effect on the bone is performed by sac- 
rificing the animal, extracting the bone therefrom and cutting the bone in order to gain ac- 
cess to the interior structure thereof. This type of method may be sufficient in order to ob- 
tain a vague indication of the effect on the bone. However, it is not desired having to 1) 

15 perform the test on animals with artificially induced conditions, 2) to sacrifice the animal 
(and thereby not having the possibility of continuing the experiment on the same animal), 
3) using a determination which requires removing the bone from the vertebrate. 

According to the invention, a technique has been found that facilitates this more desired 
20 type of investigation. Also, it is now possible to both differentiate between different bone 
defects or conditions and to determine which substance or group of substances would act 
to counteract that defect of condition. Thus, a much more specific and precise method has 
been found - and a method which may be performed fully non-invasively and without de- 
stroying the bone and without harming the vertebrate. This also provides a better manner 
25 of investigating the effects of the substances in that inter-vertebrate differences are now 
less prominent in that measurements on the same vertebrate may be performed in the 
course of the investigation. 

In a first aspect, the invention relates to a method of determining which substance, of a 
30 number of substances having an effect on bone, to administer to a vertebrate in need 
thereof, the method comprising: 

providing image data relating to one or more bones of the vertebrate, 
deriving, from the image data, one or more parameters relating to the bone, 
providing first information relating to the one or more parameters of the one or 
35 more bones of vertebrates classified as not suffering from bone diseases, 
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comparing the first information and the derived one or more parameters and de- 
termining differences there between, 

providing second information relating to one or more substances affecting bones 
and information relating to which of the one or more parameters of bones are affected 
5 thereby, 

determining which substance/substances to administer to the vertebrate on the 
basis of the result of the comparison and the second information. 

An alternative to this first aspect is one relating to an apparatus for providing information 
10 relating to the effect of a substance on a vertebrate, the apparatus comprising: 

means for providing image data relating to one or more bones of the vertebrate, 
means for deriving, from the image data, one or more parameters relating to the 

bone, 

means for providing first information relating to the one or more parameters of the 
15 one or more bones of vertebrates classified as not suffering from bone diseases, 

means for comparing the first information and the derived one or more parameters 
and determining differences there between, 

means for providing second information relating to each of the number of 
substances and information relating to which of the one or more parameters of bones are 
20 affected thereby, 

means for determining which substance/substances to administer to the vertebrate 
on the basis of the result of the comparison and the second information. 

Naturally, this apparatus may comprise means for performing any of the below mentioned 
25 steps and manners of processing the data and determining features therefrom. 

In the present context, "in need thereof would normally mean that the vertebrate has 
been diagnosed, by a physician or another expert, to not be in an optimal condition, such 
as if having a non-optimum bone mineral content, bone mineral density, fracture risk, 
30 bone mineral content, bone strength, bone quality or the like. 

Thus, preferably, the method is performed on a vertebrate or vertebrates classified as 
suffering from Bone disorders or conditions, such as type I and II osteoporosis, 
glucocroticoid induced osteoporosis and other types of drug induced osteoporosis, 
35 osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, diabetes 
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(IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone diseases such 
as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, such as breast 
cancer or prostate cancer. 

5 Conversely, bones "not suffering from bone diseases" would normally be bones, the con- 
dition (such as bone mineral density, fracture risk, bone mineral content, bone strength, 
bone quality or the like) of which is normal taking into account the weight, height, sex, 
ethnical background, age, menopausal status, life style and habits, family and/or medical 
history of the vertebrate. 

10 

Also, the first information would normally indicate normal or healthy values for the one or 
more parameters. Naturally, these values may vary between individuals and groups of 
individuals, such as dependent on weight, height, sex, ethnical background, age, meno- 
pausal status, life style and habits, family and/or medical history. It should be noted that 
15 this information might be derived from individuals, which are, in fact, not healthy, as long 
as they are not diagnosed as having bone diseases or conditions. 

"not having received the substance" will normally mean that the substance has not been 
administered in amounts where the effect, which the vertebrate is not expected to show, is 
20 not pronounced - or preferably even detectable. Preferably, these vertebrates have not 
received the substance at all. 

Normally the information relating to which of the one or more parameters are affected by 
which substances will at least comprise, for each parameter or each substance, the iden- 
25 tity of which substance is affected or which parameter it affects, respectively. However, 
below is given an example of further information, which may be provided. 

Also, the determination step may comprise deriving information on suitable doses or 
intervals of administration. 

30 

Preferably the deriving step comprises deriving, from the image data, one or more pa- 
rameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

35 - a variation of a density of cortical and/or trabecular bone of the bone. 
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The present physical distances in the bone may be any physical distance determinable 
from the image data. However, it should be noted that a distance determinable in the im- 
age data might not relate to an easily determinable distance in the bone in that the gen- 
eration of the image data may alter distances in the bone. An example of this fact is the 
5 two-dimensional image data generated by a standard X-ray image of the bone. In this im- 
age, a distance may not correlate to the same distance in the bone in that the X-ray image 
is a projection of the three-dimensional bone structure onto two dimensions. In two-di- 
mensional data, the physical distance would be a distance of e.g. a hard copy thereof or a 
metrical distance determinable between e.g. different parts of the image data. 

10 

A variation of a density of cortical and/or trabecular bone may be caused by a number of 
effects, such as the function of osteoclasts and osteoblasts in bone. These variations may 
be determinable from the image data in a number of manners. In the example where the 
image data are represented in grey values, these variations are determinable as grey 
15 value variations, such as areas having grey values different from surrounding areas. 
Again, such variations in the bone may not correspond directly to variations in the image 
data. However, a correlation will exist. 

Bone diseases/effects/conditions may exist which affect a given parameter in one direc- 
20 tion, and others may exist which affect the parameter in the opposite direction. These 
facts may point at different problems and illnesses and may be counteracted by different 
substances. Therefore, the second information may further comprise, for each substance 
and for at least one numerical parameter, information relating to whether the parameter is 
increased or decreased by the effect of the substance. In that situation, the given pa- 
25 rameter(s) would normally comprise the at least one parameter and wherein the compar- 
ing means also comprise means for deriving information relating to whether a desired ef- 
fect of the substance is to reduce or increase the parameter. 

In this context, a numerical parameter will be a parameter, which is quantifiable with e.g. a 
30 number or digit. This will normally be the case with parameters derived by computers or 
other automatic measures. 

The "desired effect" of the substance will be the effect desired in order to bring at least 
some of the one or more parameters to or toward the parameters of healthy vertebrates. 
35 In this manner, it will be determinable whether the at least one numerical parameter 
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should be decreased or increased in order to bring that or those parameter(s) to or toward 
those of healthy vertebrates. 

Thus, the determining means may comprise means for determining a combination of sub- 
5 stances providing the desired effect on the bone, where the determining means are 
adapted to take into account the information relating to the numerical parameter. 
Thus, when determining which substance or substances to administer, a single substance 
may not exist (or has not yet been identified) which counteracts the offset (difference from 
the parameter(s) of the first information) of the parameters sufficiently. Also, such a sub- 
10 stance may exist but have the "side effect" that it also has an effect on a parameter, which 
is not to be affected in the bone. In those situations, a group of substances may be identi- 
fied, the total effect of which is that desired and where e.g. "side effects" of one substance 
is counteracted by the effect of another. 

15 The determining means may be adapted to provide, in relation to at least some of the one 
or more parameters, information relating to an uncertainty, such as a statistical uncer- 
tainty and/or an uncertainty deriving from the providing of the image data or the parame- 
ters, in the deriving of the parameter(s) and to, on the basis of that uncertainty, determine 
only parameter(s) significant in relation to their related uncertainties. Even though a pa- 

20 rameter seems to vary quite significantly with a given effect, this may not be true taking 
the uncertainty into account. Thus, parameters, which are made less useful due to a large 
uncertainty, are then automatically not taken into account. 

In a first situation, the step of providing second information relating to one or more sub- 
25 stances affecting bones and information relating to which of the one or more parameters 
of bones are affected thereby comprises, for a substance: 

providing first image data relating to at least part of a bone of a vertebrate at a first 
point in time, 

administering the substance to the vertebrate after the first point in time, 
30 - providing, at a second point in time subsequent to the administering of the sub- 
stance, second image data relating to the at least part of the bone, 

deriving, from the first image data, one or more parameters relating to the bone at 
the first point in time, 

deriving, from the second image data, the one or more parameters relating to the 
35 bone at the second point in time, 
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comparing the one or more parameters relating to the first and second points in 
time, and 

providing, from the result of the comparison, information relating to which of the 
one or more parameters is affected by the substance. 

5 

Normally, the providing step will provide information relating to parameter(s) which, from 
the comparing step, is/are seen to vary with the administering of the substance. As is de- 
scribed above, this step may comprise evaluating a statistical uncertainty in order to pro- 
vide information only on parameter(s) statistically significant. 

10 

Thus, image data are provided at different points in time in order for the determination to 
be able to build on measurements performed on the same vertebrate. Preferably, the im- 
age data relate to the same bone in the vertebrate. However, it is normally contemplated 
that at least a number of bone diseases or conditions act equally on all bones in a verte- 
15 brate, whereby different image data of different bones may, in fact, provide the information 
required. 

Again, preferably, the deriving steps comprise deriving, from the image data, one or more 
parameters relating to: 

20 - one or more physical distances in the bone or the image data relating to the bone 
and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

The method of the first situation preferably further comprises administering the substance 
25 to the vertebrate a number of times between the first and second points in time and at 
predetermined doses and time intervals in order to ensure reproducible results. 

Also, in order to be able to better follow the course of action in the bone, the method may 
further comprise providing third image data of the at least part of the bone and at least 
30 one third point in time between the first and second points in time, deriving, from the third 
image data, the one or more parameters, and comparing, in the comparing step, the one 
or more parameters relating to the first, second, and third points in time. 

In order to increase the statistical certainty of the method, the image data may be pro- 
35 vided of at least parts of bones of a number of vertebrates, where the substance is pro- 
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vided to the number of vertebrates, where the one or more parameters are derived from 
each of the vertebrates, and wherein the parameters of the vertebrates are compared in 
the comparing stepin order to provide information relating to an overall effect of the 
substance on bone statistical analysis of the parameters of the vertebrates. 

5 

For a number of known parameters, only image data relating to a part of the bone suffices 
in order to derive a value for each parameter. 

A large number of bone disorders or conditions exist, and it is contemplated that some of 
10 these act differently on bones in the vertebrate and others more homogeneously. Nor- 
mally, the vertebrate or vertebrates would be classified as suffering from bone disorders 
or conditions, such as type I and II osteoporosis, glucocorticoid induced osteoporosis and 
other types of drug induced osteoporosis, osteomalacia, hyper- and hypoparathyroidism, 
hyper- and hypothyroidism, diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid 
15 arthritis, genetic bone diseases such as osteogenesis imperfecta, Paget's bone disease, 
myeloma or cancer, such as breast and/or prostate cancer. 

In a second situation, the step of providing second information relating to one or more 
substances affecting bones and information relating to which of the one or more parame- 
20 ters of bones are affected thereby comprises, for a substance: 

providing a first group of vertebrates having received the substance, 
providing a second group of vertebrates not having received the substance to any 
significant degree, 

providing a first number of image data relating to at least part of the bone of each 
25 of the first group of vertebrates, 

providing a second number of image data relating to the at least part of the bone 
of each of the second group of vertebrates, 

deriving, from the first number of image data, one or more parameters relating to 
the bones of the first group of vertebrates, 
30 - deriving, from the second image data, the one or more parameters relating to the 
bones of the second group of vertebrates, 

comparing the one or more parameters relating to the first and second groups of 
vertebrates and 

determining, from the result of the comparison, the effect of the substance on the 
35 bones of the vertebrates of the first group. 



WO 00/72216 



12 



PC17DK00/00273 



In the present context, "not having received the substance to any significant degree" will 
normally mean that the substance has not been administered in amounts where the ef- 
fect, which the vertebrate is not expected to show, is not pronounced - or preferably even 
5 detectable. Preferably, these vertebrates have not received the substance at all. 

Also in this situation is it preferred, that the deriving steps comprise deriving, from the im- 
age data, one or more parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

10 and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

Thus, instead of providing time-related information relating to e.g. a single vertebrate, in- 
formation may be provided from a number of vertebrates where some have and some 
1 5 have not received the substance - but at the same point in time. 

As described above, preferably the same bone is used in all vertebrates. However, this 
may not be required in all situations. 

20 Different types of investigations of substances exist, in one of which the vertebrates of the 
first group of vertebrates would be classified to have or have been brought to have a 
given disease or condition or a disease or condition within a given group of diseases or 
conditions, and wherein the vertebrates of the second group of vertebrates would be clas- 
sified to not have the given disease or condition or a disease or condition within the given 

25 group of diseases or conditions. 

The present classification will normally be performed by professionals, such as doctors or 
physicians. Also, professionals will typically be those actually bringing vertebrates to such 
states or conditions. It is a widely used method to bring a vertebrate to a condition being 
30 or closely resembling a given condition in order to test the effect of a given substance. 

In fact, normally, the vertebrates, typically from animals, will not actually be brought to 
have the actual disease but will be brought to have what is meant to be a model thereof. 
Thus, instead of generating the actual disease, which may, in fact, be difficult, the verte- 
35 brate is brought to a state, which is believed to be a good approximation. 



WO 00/72216 



13 



PCT/DK00/00273 



In another type of investigation, the vertebrates of the first and second groups of verte- 
brates would be classified to have a given disease or condition or a disease or condition 
within a given group of diseases or conditions - and in that situation, the method may 
comprise: 

5 - providing a third group of vertebrates, and wherein the vertebrates of the third 
group of vertebrates have been classified to not have the given disease or condition or a 
disease or condition within the given group of diseases or conditions, 

providing a third number of image data relating to the bone of each of the third 
group of vertebrates, 
10 - deriving the one or more parameters from the third image data, and 

comparing, in the comparing step, the one or more parameters of the first, second, 
and third groups of vertebrates. 

Normally, the comparing step will comprise comparing the parameters for each vertebrate 
15 individually. 

In any case, the disease or condition or group of diseases or conditions may be chosen 
from: Bone disorders or conditions, such as type I and II osteoporosis, glucocroticoid in- 
duced osteoporosis and other types of drug induced osteoporosis, osteomalacia, hyper- 
20 and hypoparathyroidism, hyper- and hypothyroidism, diabetes (IDDM), hypopituitarism, 
osteoarthritis, rheumatoid arthritis, genetic bone diseases such as osteogenesis 
imperfecta, Pagefs bone disease, myeloma or cancer, such as breast and/or prostate 
cancer. 

25 Also, in order to be able to follow the development over time, the method may comprise 
providing the first, second and optionally third sets of image data at at least two different 
points in time between which the first group has received at least part of the substance 
received, and wherein the comparing step comprises comparing the image data obtained 
at the at least two points in time. 

30 

In a third situation, the step of providing second information relating to one or more sub- 
stances affecting bones and information relating to which of the one or more parameters 
of bones are affected thereby comprises, for a substance: 

providing first image data relating to at least part of a bone of a vertebrate at a first 
35 point in time, 
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administering the substance to the vertebrate after the first point in time, 

providing, at a second point in time subsequent to the providing of the drug, sec- 
ond image data relating to the at least part of the bone, 

deriving, from the first image data, one or more parameters relating to the bone at 
5 the first point in time, 

deriving, from the second image data, the one or more parameters relating to the 
bone at the second point in time, 

comparing the one or more parameters relating to the first and second points in 

time, 

10 - determining, from the result of the comparison, whether the substance has any 
effect on the bone and 

providing, if the substance has an effect on the bone, information that the sub- 
stance has an effect on the bone and which of the one or more parameters is affected 
thereby. 

15 

Again, preferably the deriving steps comprise deriving, from the image data, one or more 
parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

20 - a variation of a density of cortical and/or trabecular bone of the bone. 

Preferably, the substance is administered to the vertebrate a number of times between 
the first and second points in time and at predetermined doses and time intervals. This is 
the most natural method of using the substance. 

25 

As is also described above, the method preferably further comprises providing third image 
data of the at least part of the bone and at at least one third point in time between the first 
and second points in time, deriving, from the third image data, the one or more parame- 
ters, and comparing, in the comparing step, the one or more parameters relating to the 
30 first, second, and third points in time. 

Also, the image data are preferably provided of at least parts of bones of a number of 
vertebrates, where the substance is provided to the number of vertebrates, where the one 
or more parameters are derived from each of the vertebrates, and wherein the parameters 
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of the vertebrates are compared in the comparing step in order to provide information 
relating to an overall effect of the substance on bone. 

In a fourth situation, the step of providing second information relating to one or more sub- 
5 stances affecting bones and information relating to which of the one or more parameters 

of bones are affected thereby comprises, for a substance: 

providing a first group of vertebrates having received the substance, 
providing a second group of vertebrates not having received the substance, 
providing a first number of image data relating to at least part of the bone of each 
1 0 vertebrate of the first group of vertebrates, 

providing a second number of image data relating to the at least part of the bone 

of each vertebrate of the second group of vertebrates, 

deriving, from the first number of image data, one or more parameters relating to 

the bones of the first group of vertebrates, 
1 5 - deriving, from the second image data, the one or more parameters relating to the 

bones of the second group of vertebrates, 

comparing the one or more parameters relating to the first and second groups of 

vertebrates, 

determining, from the result of the comparison, whether the substance has any 
20 effect on the bones of the vertebrates of the first group and 

providing, if the substance has an effect on the bone, information that the sub- 
stance has an effect on the bone and which of the one or more parameters is affected 
thereby. 

25 Again, the vertebrates of the first group of vertebrates may have been classified to have a 
given disease or condition or a disease or condition within a given group of diseases or 
conditions, and wherein the vertebrates of the second group of vertebrates have been 
classified to not have the given disease or condition or a disease or condition within the 
given group of diseases or conditions. 

30 

Alternatively, the vertebrates of the first and second groups of vertebrates may have been 
classified to have a given disease or condition or a disease or condition within a given 
group of diseases or conditions, and in that situation, the method preferably further com- 
prises: 
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providing a third group of vertebrates, and wherein the vertebrates of the third 
group of vertebrates have been classified to not have the given disease or condition or a 
disease or condition within the given group of diseases or conditions, 

providing a third number of image data relating to the bone of each of the third 
5 group of vertebrates, 

deriving, from the third image data, the one or more parameters, and 

comparing, in the comparing step, the one or more parameters of the first, second, 
and third groups of vertebrates. 

10 Also, in a second aspect, the invention relates to a method of differentiating between dif- 
ferent types of bone lesions/defects/conditions, the method comprising: 

a) providing image data on a bone in at least one vertebrate, which is suspected to suffer 
from a first bone disease or condition, 

b) providing image data on a bone in at least one vertebrate, which is suspected to suffer 
15 from a second bone disease or condition 

c) deriving one or more parameters relating to the bones based on said image data pro- 
vided in step a) and b), 

d) providing information relating to said one or more parameters of bones of vertebrates 
not suffering from the bone disease or condition, 

20 e) comparing the one or more parameters obtained in step b) with the one or more pa- 
rameters obtained in step c) to determine differences in said parameters, and 
f) identifying, on the basis of the information provided in step d) and the result of the com- 
parison of step e), one or more parameters relating to said bones of said vertebrates and 
being correlated to one or each of the first and second bone diseases or conditions. 

25 

That a vertebrate is "suspected" to suffer from a disease or condition will normally be de- 
rived by a professional, such as a doctor or a physician. The same will be the identifica- 
tion of vertebrates "not suffering" from the disease or condition. 

30 "Being correlated to" will in this respect mean that a parameter "varies with the disease or 
condition": if the disease or condition is absent, the parameter is the same for that verte- 
brate and a healthy vertebrate - if not, it is different. 
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It is clear from the above that the present invention may not only provide useful informa- 
tion relating to how to counteract certain bone diseases, it also provides useful information 
relating to exactly how substances affect bones. It will be natural to compile such informa- 
tion in order to have a much more efficient use thereof. 

5 

In a third aspect, the invention relates to a method of determining the effect of a sub- 
stance on a bone of a vertebrate, the method comprising: 

providing first image data relating to at least part of the bone at a first point in time, 

administering the substance to the vertebrate after the first point in time, 
10 - providing, at a second point in time subsequent to the administering of the sub- 
stance, second image data relating to the at least part of the bone, 

deriving, from the first image data, one or more parameters relating to the bone at 
the first point in time, 

deriving, from the second image data, the one or more parameters relating to the 
15 bone at the second point in time, 

comparing the one or more parameters relating to the first and second points in 
time, and 

determining, from the result of the comparison, the effect of the substance on the 

bone 

20 

wherein the deriving steps comprise deriving, from the image data, one or more parame- 
ters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

25 - a variation of a density of cortical and/or trabecular bone of the bone. 

An alternative to this third aspect is one relating to an apparatus for determining the effect 
of a substance on a bone of a vertebrate, the apparatus comprising: 

means for providing first image data relating to at least part of the bone at a first 
30 point in time, 

means for providing, at a second point in time, second image data relating to the 
at least part of the bone, 

means for deriving, from the first image data, one or more parameters relating to 
the bone at the first point in time, 
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means for deriving, from the second image data, the one or more parameters 
relating to the bone at the second point in time, 

the deriving means being adapted to derive, from the image data, one or more 
5 parameters relating to: 

one or more physical distances in the bone or the image data relating to 
the bone and/or 

a variation of a density of cortical and/or trabecular bone of the bone 
means for comparing the one or more parameters relating to the first and second 
10 points in time, and 

means for determining, from the result of the comparison, the effect of the 
substance on the bone. 

Preferably, naturally, the first image data is provided at a point in time prior to that of the 
15 providing the second image data - and preferably, the bone has received the substance at 
at least one point in time between these two points in time. 

This apparatus may comprise means for performing any of the below method steps 
related to the processing of the data. 

20 

Thus, image data are provided at different points in time in order for the determination to 
be able to build on measurements performed on the same vertebrate. Preferably, the im- 
age data relate to the same bone in the vertebrate. However, it is normally contemplated 
that at least a number of bone diseases or conditions act equally on all bones in a verte- 
25 brate, whereby different image data of different bones may, in fact, provide the information 
required. 

The method preferably further comprises administering the substance to the vertebrate a 
number of times between the first and second points in time and at predetermined doses 
30 and time intervals in order to ensure reproducible results. 

Also, in order to be able to better follow the course of action in the bone, the method may 
further comprise providing third image data of the at least part of the bone and at least 
one third point in time between the first and second points in time, deriving, from the third 
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image data, the one or more parameters, and comparing, in the comparing step, the one 
or more parameters relating to the first, second, and third points in time. 

In order to increase the statistical certainty of the method, the image data may be pro- 
5 vided of at least parts of bones of a number of vertebrates, where the substance is pro- 
vided to the number of vertebrates, where the one or more parameters are derived from 
each of the vertebrates, and wherein the parameters of the vertebrates are compared in 
the comparing step. 

10 For a number of known parameters, only image data relating to a part of the bone suffices 
in order to derive a value for each parameter. 

A large number of bone disorders or conditions exist, and it is contemplated that some of 
these act differently on bones in the vertebrate and others more homogeneously. Nor- 

15 mally, the vertebrate or vertebrates would be classified as suffering from bone disorders 
or conditions, such as type I and II osteoporosis, glucocorticoid induced osteoporosis and 
other types of drug induced osteoporosis, osteomalacia, hyper- and hypoparathyroidism, 
hyper- and hypothyroidism, diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid 
arthritis, genetic bone diseases such as osteogenesis imperfecta, Paget's bone disease, 

20 myeloma or cancer, such as breast and/or prostate cancer. 

In yet another aspect, the invention relates to a method of determining the effect of a sub- 
stance on a bone of a vertebrate, the method providing: 

providing a first group of vertebrates having received the substance, 
25 - providing a second group of vertebrates not having received the substance to any 
significant degree, 

providing a first number of image data relating to at least part of the bone of each 
of the first group of vertebrates, 

providing a second number of image data relating to the at least part of the bone 
30 of each of the second group of vertebrates, 

deriving, from the first number of image data, one or more parameters relating to 
the bones of the first group of vertebrates, 

deriving, from the second image data, the one or more parameters relating to the 
bones of the second group of vertebrates, 
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comparing the one or more parameters relating to the first and second groups of 
vertebrates and 

determining, from the result of the comparison, the effect of the substance on the 
bones of the vertebrates of the first group 

5 

wherein the deriving steps comprise deriving, from the image data, one or more parame- 
ters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

10 - a variation of a density of cortical and/or trabecular bone of the bone. 

Thus, instead of providing time-related information relating to e.g. a single vertebrate, in- 
formation may be provided from a number of vertebrates where some have and some 
have not received the substance - but at the same point in time. 

15 

In the present context, "not having received the substance" will mean that the substance 
has not been administered in amounts where the effect, which the vertebrate is not ex- 
pected to show, is pronounced - or even detectable. Preferably, these vertebrates have 
not received the substance at all. 

20 

As described above, preferably the same bone is used in all vertebrates. However, this 
may not be required in all situations. 

Different types of investigations of substances exist, in one of which the vertebrates of the 
25 first group of vertebrates would be classified to or have been brought to have a given dis- 
ease or condition or a disease or condition within a given group of diseases or conditions, 
and wherein the vertebrates of the second group of vertebrates would be classified to not 
have the given disease or condition or a disease or condition within the given group of 
diseases or conditions. 

30 

In fact, normally, the vertebrates, typically from animals, will not actually be brought to 
have the actual disease but will be brought to have what is meant to be a model thereof. 
Thus, instead of generating the actual disease, which may, in fact, be difficult, the verte- 
brate is brought to a state, which is believed to be a good approximation. 

35 
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In another type of investigation, the vertebrates of the first and second groups of verte- 
brates would be classified to have a given disease or condition or a disease or condition 
within a given group of diseases or conditions - and in that situation, the method may 
comprise: 

5 - providing a third group of vertebrates, and wherein the vertebrates of the third 
group of vertebrates have been classified to not have the given disease or condition or a 
disease or condition within a given group of diseases or conditions, 

providing a third number of image data relating to the bone of each of the third 
group of vertebrates, 
10 - deriving the one or more parameters from the third image data, and 

comparing, in the comparing step, the one or more parameters of the first, second, 
and third groups of vertebrates. 

In any case, the disease or condition or group of diseases or conditions may be chosen 
15 from: Bone disorders or conditions, such as type I and II osteoporosis, glucocroticoid in- 
duced osteoporosis and other types of drug induced osteoporosis, osteomalacia, hyper- 
and hypoparathyroidism, hyper- and hypothyroidism, diabetes (IDDM), hypopituitarism, 
osteoarthritis, rheumatoid arthritis, genetic bone diseases such as osteogenesis 
imperfecta, Paget's bone disease, myeloma or cancer, such as breast and/or prostate 
20 cancer. 

Also, in order to be able to follow the development over time, the method may comprise 
providing the first, second and optionally third sets of image data at at least two different 
points in time between which the first group has received at least part of the substance 
25 received, and wherein the comparing step comprises comparing the image data obtained 
at the at least two points in time. 

Also, in one aspect, the invention relates to a type of Drug Discovery, that is, a method of 
screening test substances to identify substances which have an effect on a bone of a 
30 vertebrate, the method comprising: 

providing first image data relating to at least part of the bone at a first point in time, 
administering the substance to the vertebrate after the first point in time, 
providing, at a second point in time subsequent to the providing of the drug, sec- 
ond image data relating to the at least part of the bone, 
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deriving, from the first image data, one or more parameters relating to the bone at 
the first point in time, 

deriving, from the second image data, the one or more parameters relating to the 
bone at the second point in time, 
5 - comparing the one or more parameters relating to the first and second points in 
time, and 

determining, from the result of the comparison, whether the substance has any 
effect on the bone 

10 wherein the deriving steps comprise deriving, from the image data, one or more parame- 
ters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

a variation of a density of cortical and/or trabecular bone of the bone 

15 

Preferably, the substance is administered to the vertebrate a number of times between 
the first and second points in time and at predetermined doses and time intervals. This is 
the most natural method of using the substance. 

20 As is also described above, the method preferably further comprises providing third image 
data of the at least part of the bone and at at least one third point in time between the first 
and second points in time, deriving, from the third image data, the one or more parame- 
ters, and comparing, in the comparing step, the one or more parameters relating to the 
first, second, and third points in time. 

25 

Also, the image data are preferably provided of at least parts of bones of a number of 
vertebrates, where the substance is provided to the number of vertebrates, where the one 
or more parameters are derived from each of the vertebrates, and wherein the parameters 
of the vertebrates are compared in the comparing step. 

30 

In a further aspect, the invention relates to a method of screening test substances to iden- 
tify substances which have an effect on a bone of a vertebrate, the method comprising: 
providing a first group of vertebrates having received the substance, 
providing a second group of vertebrates not having received the substance, 
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providing a first number of image data relating to at least part of the bone of each 
vertebrate of the first group of vertebrates, 

providing a second number of image data relating to the at least part of the bone 
of each vertebrate of the second group of vertebrates, 
5 - deriving, from the first number of image data, one or more parameters relating to 
the bones of the first group of vertebrates, 

deriving, from the second image data, the one or more parameters relating to the 
bones of the second group of vertebrates, 

comparing the one or more parameters relating to the first and second groups of 
10 vertebrates and 

determining, from the result of the comparison, whether the substance has any 
effect on the bones of the vertebrates of the first group 

wherein the deriving steps comprise deriving, from the image data, one or more parame- 
15 ters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

20 Again, the vertebrates of the first group of vertebrates may have been classified to have a 
given disease or condition or a disease or condition within a given group of diseases or 
conditions, and wherein the vertebrates of the second group of vertebrates have been 
classified to not have the given disease or condition or a disease or condition within the 
given group of diseases or conditions. 

25 

Alternatively, the vertebrates of the first and second groups of vertebrates may have been 
classified to have a given disease or condition or a disease or condition within a given 
group of diseases or conditions, and in that situation, the method preferably further com- 
prises: 

30 - providing a third group of vertebrates, and wherein the vertebrates of the third 
group of vertebrates have been classified to not have the given disease or condition or a 
disease or condition within the given group of diseases or conditions, 

providing a third number of image data relating to the bone of each of the third 
group of vertebrates, 

35 - deriving, from the third image data, the one or more parameters, and 
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comparing, in the comparing step, the one or more parameters of the first, second, 
and third groups of vertebrates. 

It may be desired that each vertebrate be a laboratory or test animal, and wherein the 
5 method comprises finally sacrificing the animal in order to e.g. be able to perform more 
advanced studies of the vertebrate subsequent to the other parts of the method. 

Independently on which of the aspects of the invention is used, a number of methods are 
known for providing image data relating to bones. Each of these methods has advantages 
10 and disadvantages - such as the obtainable resolution, the minimum radiation dose re- 
quired, etc. 

It is presently preferred to acquire information relating to the internal structure of the bone 
- and in a manner where the bone is not destroyed, the image data are preferably ob- 
15 tained by a non-invasive method, such as by exposing the bone or bones to radiation, 
such as X-rays, performing a CT scanning, a NMR scanning or exposing the bone or 
bones to ultra sound. 

In order to obtain reasonable calculation speeds, the image data are preferably two-di- 
20 mensional image data relating to at least part of the bone or bones. Also, depending on 
the actual method, only a part of a single bone, such as ultra sound on the heel, is used, 
and in other methods, such as X-ray of a hand, several bones may be seen in the image 
data. 

25 A vast number of parameters may be derived from the above methods. Parameters of 
these types may be seen in e.g. WO 96/07161 , WO 97/42602, and WO 99/01835 - all of 
which are hereby incorporated by reference. 

One situation is one where at least one of the one or more parameters relates to a physi- 
30 cal distance within the bone. It is well known to derive such parameters from e.g. X-ray 
images. 
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In that situation, the image data may comprise information relating to the cortical bone of 
at least a part of the bone, and one of the one or more parameters may relate to a thick- 
ness, t, of the cortical bone. In addition or alternatively, one of the one or more parame- 
ters may relate to a width, w, of the cortical bone. 

5 

These thicknesses may be derived from e.g. metacarpal bones - and even from different - 
or both - sides thereof. However, this is not desired with all bones. If the bone is a radius, 
the t value or values is/are preferably determined on a radial side of the radius, and if the 
bone is an ulna, the t value or values is/are preferably determined on an ulnar side of the 
10 ulna. 

t and w may be determined along a single line extending in a direction perpendicular to a 
longitudinal axis of the bone, and preferably pairs of (t,w) are determined for a plurality of 
lines extending in the direction and being positioned at different positions along the longi- 
15 tudinal direction of the bone, and at least one of the one or more parameters may be de- 
termined on the basis of pairs of (t,w) corresponding to the individual lines. In that situa- 
tion, the at least one parameter may be determined on the basis of mean values of the t 
and w values corresponding to the individual lines. 

20 In order to provide a better measurement, t and w values may be determined for each of 
more than 10, such as more than 20, such as more than 40, such as more than 60, such 
as more than 80, such as more than 1 00, such as more than 110 lines per cm. of the 
bone within a predetermined longitudinal part of the bone. 

25 In another situation, at least one of the one or more parameters is derived on the basis of 
a determination of at least one value representing a variation of a density of cortical 
and/or trabecular bone of the at least part of the bone. 

For example, the at least one variation value may be related to a cortical porosity of the 
30 bone. 

Then, the at least one variation value may be determined by a method comprising deter- 
mining, within the image data, first parts thereof relating to cortical bone and second parts 
positioned within the first parts, the second parts not relating to cortical bone. The pa- 
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rameter may then relate to e.g. a relationship between those parts, such as the relative 
sizes or positioning thereof. 

When the image data are derived from a digitised, two-dimensional image of the at least 
5 part of the bone, the at least one variation value may be derived by determining one or 
more first areas in the image data relating to or representing cortical tissue and one or 
more second areas positioned within the first area(s), each of the second area(s) repre- 
senting a local deviation in grey value. Then the local deviation in grey value may be a 
grey value being lower than the grey values of any surrounding areas of the first area(s), 
1 0 such as lower than the grey values of the surrounding areas minus a predetermined 
value. 

When the at least one variation value is related to a cortical porosity or a variation in the 
cortical part of the image data, the at least one variation value may be determined by, in 
15 the image data: 

determining an outer cortical edge of the bone and denoting pixels thereon as cor- 
tical pixels, 

repeating: 

for each, current cortical pixel, denote any of its neighbouring pixels having 
20 a grey value higher than that of itself as a cortical pixel, 

until no non-cortical pixels in the neighbourhood of any current cortical 
pixel have a grey level value that is greater than the grey level value of the given current 
cortical pixel, thereby defining an inner cortical edge as the inner (endosteal) boundary of 
the final set of cortical pixels, 
25 - denoting non-cortical image data pixels within the outer and inner cortical edges 
as porosity pixels, and 

deriving the at least one variation value on the basis of the cortical pixels and the 
porosity pixels. 

30 Thus, a procedure is run where a single pixel on the boundary of the bone will "send out 
scouts" that also send out scouts in directions where the grey values increase. In that 
manner, areas with low grey values are identified. 
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Thus, the deriving step may comprise deriving the at least one variation value on the ba- 
sis on a number of porosity pixels compared to a number of cortical pixels in the image 
data, and the deriving step may comprise deriving the at least one variation value on the 
basis of a number, size, and shape of groups of porosity pixels. Also, the groups of poros- 
5 ity pixels may be isolated groups of porosity pixels being surrounded by non-porosity pix- 
els. 

When the at least one variation value is related to a cortical porosity or a variation in the 
cortical part of the image data, the at least one variation value may alternatively or addi- 
10 tionally be determined from oblong density variations in the image data, such as of parts 
of the image data relating to cortical tissue, such as where the at least one value is de- 
termined from oblong density variations extending at least substantially in a longitudinal 
direction of the bone. 

15 These oblong variations are seen in both healthy and osteoporotic bone, and preferably, 
the oblong density variations have a length to width ratio of at least 1%:1 , such as at least 
2:1, such as at least 2!4:1, preferably at least 3:1, such as at least 3.5:1, preferably at 
least 4:1, such as at least 4.5:1. However, for osteoporotic bones, it may be desired to 
actually have these longer, as investigations point to the fact that these channels in bone 

20 tend to be longer in osteoporotic bones. 

The oblong density variations may have a length to width ratio in the interval 1 1 >4-8:1, such 
as 1.7-6:1, preferably 2-4:1. Also, it may be preferred to increase the length to width ratio 
due to the above-described difference between the channels in healthy and osteoporotic 
25 bones. 

Normally, the oblong density variations in the image data will relate to oblong density 
variations, such as cavities, tunnels, channels, valleys, grooves, or the like, in the cortical 
bone. 

30 

Normally, the density variations of the bone have a length in the interval of 0.4-15 mm, 
such as in the interval 0.5-12 mm, preferably in the interval 0.8-10 mm, such as in the in- 
terval 1-8 mm, such as on the order of 2 mm, and a width in the interval 0.1-6 mm, such 
as in the interval 0.13-5 mm, preferably in the interval 0.15-3 mm, such as in the interval 
35 0.17-2 mm, preferably in the interval 0.18-1 mm, such as on the order of 0.2-0.4 mm. 
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However, a longer length may be expected in bones diagnosed with certain types of dis- 
eases or conditions. 

Also, the oblong density variations of image data may have a length in the interval of 0.4- 
5 15 mm, such as in the interval 0.5-12 mm, preferably in the interval 0.8-10 mm, such as in 
the interval 1-8 mm, such as on the order of 2 mm, and a width in the interval 0.1-6 mm, 
such as in the interval 0.13-5 mm, preferably in the interval 0.15-3 mm, such as in the in- 
terval 0.17-2 mm, preferably in the interval 0.18-1 mm, such as on the order of 0.2-0.4 
mm. Naturally, this depends on the magnification in the image generating set-up. Also, 
10 slightly longer lengths may be suitable at least for bones with certain types of osteoporo- 
sis. 

In one situation, the determined variation value is determined from a number of oblong 
density variations in the image data, a mean width thereof, and/or a mean length thereof. 

15 

Additionally or alternatively, the determined variation value may be determined from a 
grey value difference between the oblong density variations and other parts of the image 
data. 

20 The variation value may also or alternatively be determined by: 
obtaining the power spectrum relating to the image data, 
identifying parts thereof relating to frequencies or frequency intervals correspond- 
ing to predetermined dimensions or dimension intervals of the oblong variations, and 

a) estimating the energy at or within the frequencies or frequency in- 

25 tervals, and/or 

b) removing at least substantially all other parts of the power spec- 
trum, inversely Fourier transforming the resulting power spectrum, and determining a 
variation of the grey levels of the resulting image. 

30 Additionally or alternatively, the variation value may be determined by filtering the image 
data with an oblong Gaussian Kernel and subsequently determining a variation of grey 
levels of the resulting image. In that situation, the variation value is preferably determined 
by filtering each two copies of the image data with an oblong Gaussian Kernel, the two 
oblong Gaussian kernels having different dimensions, and subsequently subtracting the 

35 two resulting image data and determining a variation of grey levels of the resulting image. 
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The oblong kernel(s) may have a length to width ratio of at least 1 14:1 , such as at least 
2:1, such as at least 214:1, preferably at least 3:1 , such as at least 3.5:1, preferably at 
least 4:1 , such as at least 4.5:1 , and a length to width ratio in the interval 114-8:1 , such as 
1 .7-6: 1 , preferably 2-4:1 . 

5 

Also, the or one oblong kernel may have a length in the interval of 0.4-15 mm, such as in 
the interval 0.5-12 mm, preferably in the interval 0.8-10 mm, such as in the interval 1-8 
mm, such as on the order of 2 mm, and a width in the interval 0.1-6 mm, such as in the 
interval 0.13-5 mm, preferably in the interval 0.15-3 mm, such as in the interval 0.17-2 
10 mm, preferably in the interval 0.18-1 mm, such as on the order of 0.2-0.4 mm - again, a 
longer kernel may be chosen for certain types of osteoporosis. 

The oblong kernel may have a length in the interval of 1-15 mm, such as in the interval 
1.3-10 mm, preferably in the interval 114-6 mm, such as in the interval 2-4 mm, such as on 
15 the order of 3 mm, and a width in the interval 0.4-6 mm, such as in the interval 34-514 mm, 
preferably in the interval 0.7-5 mm, such as in the interval 1-2!4 mm, preferably in the in- 
terval 1 .2-2 mm, such as on the order of 114 mm - or slightly longer. 

Also, the at least one variation value may determined by: 
20 - providing a pre-determined template of the oblong density variations, 

template matching parts of the image data with the predetermined shape so as to 
identify oblong density variations, 

determining the variation value from: 

the number of identified oblong density variations, 
25 - a grey value difference between the identified oblong density variations 

and other parts of the image data, 

an average length and/or width of the identified oblong density variations. 

In fact, as is also described above, different types of osteoporosis or bone disorders or 
30 defects in general are expected to generate different types of alterations of the internal 
structure of bones. Therefore, such a difference is contemplated to be a change of the 
dimensions of the oblong density variations of the cortical part of the image data. Thus, by 
looking for different dimensions, information may be derived as to which type of disorder 
or defect has actually generated the problem - and also as to which substance(s) coun- 
35 teract the actual problem. Thus, parameters may relate to different dimensions, or the im- 
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age data may be scanned with different kernel dimensions in order to determine the best 
fit - and thereby the actual dimensions. 

In another situation, the image data comprises information relating to the trabecular 

5 structure of at least a part of a bone of the vertebrate, and the at least one variation value 
is derived by one or more of the following methods: 

obtaining an estimate of the parametric estimate of the power spectrum of the im- 
age data and extracting at least one of the one or more parameters as one relating to the 
energy distribution of the parametric estimate, 

10 - obtaining, on the basis of image data on which a Fourier method has been used to 
emphasise the information in the image data relating to the trabecular structure, an esti- 
mate of a grey-level co-occurrence matrix and extracting at least one of the one or more 
parameters on the basis of the estimated co-occurrence matrix, 

obtaining an estimate of the projected trabecular pattern of the image data by us- 

15 ing a Fourier method to emphasise the information in the image data relating to the trab- 
ecular structure and subjecting the manipulated image data to a morphological operation, 
and extracting at least one of the one or more parameters as one relating to the trabecular 
structure from the estimated projected trabecular pattern, and 

obtaining, on the basis of a frequency analysis of the image data, at least one of 

20 the one or more parameters as one relating to the periodicity of the trabecular structure of 
the part of the bone. 

A vast number of different parameters may be generated in this manner. These are de- 
scribed in detail in WO 96/07161 . 

25 

Also, the image data may comprise information relating to the trabecular structure of at 
least a part of a bone of the vertebrate, and wherein the at least one variation value may 
be derived by: 

determining values reflecting the projected trabecular density in the image data, 
30 caused by the X-ray attenuating properties of cancellous bone in the part of the bone, for 
each of a number of locations or areas in the image data, and 

deriving the at least one value from the variation of the determined PTD-values, 
preferably in the longitudinal direction of the bone. 
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In this manner, a density profile may be derived where from a number of parameters may 
be derived. This is extensively described in WO 97/42602. 

Normally, when parameters are compared relating to different vertebrates or bones with 
5 different conditions, the method may include determining, from the result of the compari- 
son, which of the one or more parameters is/are effected by the substance. 

Also, it is preferred that a statistical uncertainty is derived in relation to each of the one or 
more parameters, and wherein the determined parameters are that/those having an effect 
1 0 exceeding the uncertainty(ies) relating to the parameters) 

The effect of the substance on the bone may relate to e.g. the change in the one or more 
parameters or overall parameters of the bone, such as BMD, strength, bending strength, 
thickness, width, weight, and/or it may relate to the changes in e.g. bone formation 

15 brought about by the substance, whereby the substance may be bone forming, anti 
resorptive or have a heterogeneous effect. The overall effect/operation of the substance 
may be one relating to the function of the osteoblasts, osteocytes, osteoclasts or other 
parts of the vertebrae or body, such as the operation of the endocrine, paracrine or 
autocrine systems, such as growth hormone or growth hormone releasing substances or 

20 hormones) Nevertheless, if this operation or function has an impact on bone, this will be 
determinable from the image data. 

Preferably, the substance is a substance known to or expected to have an effect on 
bones. Substances of this type are: 
25 - an antiresorptive agent, such as an agent expected to act through a direct or indi- 
rect inhibition of the osteoclasts, such as: 

Bisphosphonates, such as Etidronate, Clodronate, Pamidronate, 
Alendronate, Risedronate, Ibandronate, Zoledronate, Tiludronate, Incadronate, 
Neridronate, Olpadronate, EB-1053, or YH 529, 
30 - Hormone Replacement Therapy , such as Estrogens, Estrogens with 

progestogens, Tibolone, or Trimegestone, 

Selective Estrogen Receptor Modulator (SERM), such as Raloxifene, 
Droloxifene, Idoxifene, Levormeloxifene, or Tamoxifen, or 

Calcitonin, such as Salmon calcitonin, Human calcitonin, Calcitonin gene- 
35 related peptide. 
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an agent with Heterogeneous Effect, such as Stanozolol, Oxandrolone, Nandro- 
lone, Dihydrotachysterol and other vitamin D derivatives, Alfacalcidiol, Calcifero, Calcitriol, 
Cholecalciferol, ST-630, Calcium, Calcium with vitamin D t Ipriflavone and other isofla- 
vones, or Strontium, 

5 - a bone forming agent, such as an agent expected to act through a direct or indirect 
stimulation of the osteoblasts through a mitogenic mechanism, such as Sodium fluoride, 
Monofluorophosphate, Parathyroid hormone^, Parathyroid hormone^, Growth hor- 
mone, Growth hormone releasing compounds, IGF-1, IGF-BP3, Osteogenic potein-1, or 
Bone morphogenetic protein-2. 

10 

Even though virtually any vertebrate may be used, the fact that the present methods fa- 
cilitate non-invasive measurement, at least the determination of which substance to ad- 
minister will normally be performed on humans, as osteoporosis is rarely treated in ani- 
mals. However, normally, the vertebrate may be that of a human, a horse, a large ape, a 
15 great ape, an anthropoid ape, a monkey, a pig, a cow, a rat, a rabbit, a dog, a cat 

Also, as the relevant parameters may relate to both cortical and trabecular tissue, virtually 
any bone in the vertebrate may be used. However, normally, the bone is taken from the 
group consisting of radius, femur, corpus vertebrae (L1, L2, L3, L4, L5, T1, T2, T3, T4, T5, 
20 T6, T7, T8, T9, T10, T11, T12, C1, C2, C3, C4, C5, C6, C7), calcaneus, talus, os carpi, 
metatars, metacarpi, phalanges, tibia, fibula, patella, ulna, humerus, mandible, clavicula, 
scapula, os coxae, os naviculare, os cuboideum, os cuneiform I, os cuneiform II, os cu- 
neiform III, os sacrum, os coccygis. 

25 One of the advantages of using especially image information relating to e.g. not only a 
single plane of the bone, as is the case when a bone has been removed, cut and the ex- 
posed structure is viewed, but rather to two-dimensional image data generated by radia- 
tion transmitted through the bone in order to obtain information relating to the structure 
throughout the thickness of the bone. In that situation, the image data will comprise more 

30 information, and the statistical uncertainty always incurred in measurements on organic 
tissue, will be reduced due to the super-positioning of the features of the bone. Also, this 
is more easily performed non-invasively. This also has the advantage that the bone need 
not be removed and destroyed in order to obtain the data. This is an advantage both seen 
from the patient, due to the time and resources required in order to remove and cut the 
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bone, and due to the fact that the same bone may be measured several times during a 
test. 

Specifically, it has been found that a parameter such as especially cortical porosity is af- 
5 fected by Bisphosphonates, and especially cortical thickness is affected by HRT. 

The invention also relates to a system for identifying substances for treatment of bone 
diseases/conditions, the system comprising: 

a data storage means comprising identifiers for a number of substances and, for 
10 each substance, data relating to parameters affected by the substance when administered 
to a vertebrate, 

search means for comparing one or more given parameter(s) to parameters in the 
data storage means in order to identify one or more substances affecting one or more 
given parameters), 
15 - means for outputting the identity of the identified substances. 

Normally, this system will be computer-based. 

Thus, when presenting the system to one or more parameters, which are desired affected, 
20 the system responds by providing at least the identity of substances found to affect such 
parameters). 

Preferably, the storage means are adapted to hold data relating to one or more bone pa- 
rameters for each substance, each parameter relating to: 
25 - one or more physical distances in the bone or the image data relating to the bone 
and/or 

a variation of a density of cortical and/or trabecular bone of the bone, and 
wherein each of the one or more given parameters is one of the plurality of parameters. 

30 The search means may be adapted to identify a combination of substances affecting the 
one or more given parameters of the bone. 

Also, the data storage is adapted to hold data relating to at least one numerical parameter 
for each substance as well as information relating to whether the at least one numerical 
35 parameter is increased or decreased by the effect of the substance. 



WO 00/72216 



34 



PCT/DKOO/00273 



In that situation, the given parameter(s) comprise(s) at least one of the at least one 
numerical parameter, and wherein the search means further comprise means for 
receiving information relating to whether a desired effect of the substance is to reduce or 
increase the at least one numerical parameter. 

5 

Thus, the search means may comprise means for determining a combination of sub- 
stances providing a desired effect on the bone, where the search means are adapted to 
take into account the information relating to the numerical parameter. This also has the 
advantage that, if such information is provided for a number of or all of the substances, to 
1 0 actually provide a specific combination of substances which take the actual parameters of 
an individual vertebrate into account and act to optimise the parameters of that individual. 
Often, the effects of a numerical parameter of a combinations of substances may simply 
be added or subtracted, taking into account different doses etc. 

15 Also, and especially when a number of substances are identified, the data storage pref- 
erably further comprises, for each substance, information relating to recommendable 
doses to be administered. Also, the data storage may further comprise, for each sub- 
stance, information relating to recommendable time periods between doses - and infor- 
mation relating to any known side effects of the substance. 

20 

In the following, a preferred embodiment of the invention is described in four examples 
and with reference to the drawing, wherein 

- a test has been performed for certain substances, and wherein 

25 - the annual changes, expressed in percentages, for the untreated group is shown in fig- 
ure 1, 

- the SD normalised annual changes for the untreated group is shown in figure 2, 

30 - the annual changes, expressed in percentages, for the "HRT > 90%" group is shown in 
figure 3, 

- the SD normalised annual changes for the "HRT > 90%" group is shown in figure 4, 
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- the annual changes, expressed in percentages, for bisphosphonate group is shown in 
figure 5, 

-the SD normalised annual changes for the bisphosphonate group is shown in figure 6, 

5 

- the actual annual changes for porosity and striation for all groups are shown in figure 7, 

-the annual changes, expressed in percentages, for the untreated and the bisphospho- 
nate treated group are shown in figure 8, 

10 

- the actual annual changes for porosity and striation for the untreated and the bisphos- 
phonate treated groups are shown in figure 9, 

-the annual changes, expressed in percentages, for the untreated and the HRT treated 
15 group are shown in figure 10, 

- the actual annual changes for porosity and striation for the untreated and the HRT 
treated group are shown in figure 11, and where 

20 - figure 12 illustrates the result from a left hand of the porosity algorithm for meta- 
carpal 4,3,2 and for ulna and radius, 

figure 1 3, illustrates the corresponding original images relating to figure 12, 

25 - figure 14 illustrates different types of bone resorption in a cross section of a radius, 

figure 15 illustrates, in a one-dimensional case, an original profile shown together 
with this profile filtered with a medium frequency filter and a low frequency filter, 

30 - figure 16 illustrates the cortical thickness t and the bone width W in a cross sec- 
tion of a cortical bone and the thickness profile of the bone when projected onto the film. 

figure 17 illustrates the regions used in a Hologic QDR 2000 DEXA apparatus as 
made by Hologic, USA, 
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figure 18 illustrates the regions of interest used in the Pronosco X-posure Sys- 
tem™ embodying the technology, 

figure 19 illustrates a system for incorporating a preferred method according to the 
5 invention. 

Example 1: determination of cortical porosity in X-ray images 

The actual generation of an X-ray image is well known to the skilled person. Digitising this 
10 image is also a well-known procedure. However, both procedures are described in the 
Applicants co-pending applications, such as WO 96/07161 , WO 97/42602, and WO 
99/01835. 

The present method derives a value for the cortical porosity from an X-ray image com- 
15 prising at least some of the metacarpals and the distal radius of a person. This part of the 
person makes it possible to use a low radiation dose and thereby facilitate screening of 
persons. 

1.1 The algorithm used 

20 

The input to the algorithm is a grey scale image of a mainly cortical bone where the peri- 
osteal boundary of the bone has already been defined. The output is a segmentation of 
the bone into cortex pixels and non-cortex pixels, i.e. each pixel in the bone is classified 
as belonging either to the cortex or not. The cortex pixels of one half of the bone are a 
25 connected region. This segmentation takes place as follows: 
Rotate the bone so that it is vertical in the image 
Define all pixels on the periosteal boundary as cortex pixels 
A pixel has eight neighbours. For every cortex pixel c, investigate these eight 
neighbours in turn (it is not allowed to go beyond the periosteal boundary, and it is not 
30 necessary to investigate pixels that are already identified as cortex pixels. If the pixel 
value of the neighbour is higher than or equal to the pixel value in pixel c, add the neigh- 
bour pixel to the list of cortex pixels. 

Repeat step 3 until no further cortex pixels can be constructed. 
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This algorithm is called the greedy algorithm due to the appetite of step 4. The porosity is 

now defined as follows: 

At every horizontal line of the image, find the midpoint of the bone 
Move towards the cortex until the first cortex pixel is met. This defines the cortical 
5 boundary. (This is done in each side, however, in the radius and ulna only in one side - 

the radial and ulnar side, respectively). The distance to the periosteal surface is called the 

cortical thickness. 

Count the number of non-cortex pixels in the range of the cortical thickness. These 
are the porosity pixels. 

10 - Compute the two-dimensional porosity p (small p) as the area of the porosity pix- 
els divided by the area defined by the cortical boundary and the periosteal boundary. This 
is done in the bones: metacarpal 2, 3, 4, ulnar side of ulna, radial side of radius. 
Define the combination 2-d porosity as 
p(combined) = 1 /* (p(radius) + 1/3 (p(met1) + p(met2) + p(met3)) 

15 

It may be desired to also have ulna take part of this calculation. 
Transform from 2-d to 3-d according to 
P = p(combined)/sqrt(T71 mm) - 2 
where 7 is the mean metacarpal thickness 
20 - Scale and shift P arbitrarily to lie within the range 1-9, and ensure that all future 
porosities are in this range by truncation to this range. 

Preferably, the present calculation is performed on the "raw" image without any un- 
sharpening or filtering - and even without a background correction removing the effect of 
25 soft tissue. 

In Figure 12 is shown the result from a left hand of the porosity algorithm for metacarpal 
4,3,2 and for ulna and radius. The black pixels are the porosities. Notice that endosteal 
"bays', which are shielded from horizontal "wind", count as porosity. This is a matter of 
30 definition and could be changed in later variants of the algorithm. 



In Figure 13, the corresponding original images are illustrated. 
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Example 2: determination of cortical striation in X-ray images 

The present methodology is preferably performed on X-ray images as those of Example 
1. 

5 

In Figure 14, different types of bone resorption are illustrated in a cross section of a ra- 
dius. 

A preferred manner of deriving the variation value and especially for deriving information 
10 relating to spongiosation or intra-cortical resorption - cortical striation - is to derive infor- 
mation relating to oblong density variations in the image data. 

Two manners are presently preferred in that respect, one of which uses a Difference Of 
Gaussians (DOG) procedure. 

15 

The Difference Of Gaussians may be seen as an equivalent to a band pass filter. In this 
procedure, two copies of the image data are filtered with two different filters - or Gaussian 
kernels having different sizes - one for each image. Naturally, this filtering will alter the 
image data. 

20 

In Figure 15, this is seen in a one-dimensional case, where an original profile is shown 
together with this profile filtered with a medium frequency filter and a low frequency filter. 

Subsequently, the two images are combined, whereby differences there between are en- 
25 hanced - differences relating to the frequency band between the two filters. 

In two dimensions, this procedure may be used with different frequencies in the two di- 
mensions. In this manner, oblong density variations may be identified, as different fre- 
quencies may be identified in the two directions. 

30 

Consequently, two copies of the image data are generated, one being filtered with a 
Gaussian kernel having a size of 4x1 mm and one having a size of 12x3 mm, the longer 
dimension being along the longitudinal direction of the bone - the predominant direction of 
the spongiosation. 

35 
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When subtracting these images, information will remain as to oblong density variations 
having dimensions within the intervals identified by the different kernels. 
Subsequent to that, the spongiosation may be quantified by counting the number of such 
density variations, quantifying a mean grey level of the image data - over a background 
5 level - (due to the fact that only these oblong density variations are represented in the im- 
age data), or quantifying dimensions thereof. 

This quantification of the density variation may be used as the variation value used in the 
present method. 

10 

Another method of deriving the information relating to these oblong density variations may 
be to Fourier transform the image data to obtain the power spectrum thereof. The image 
data being two dimensional, the power spectrum is also two dimensional. Therefore, 
points or areas may be defined at non-equal frequencies for the two dimensions, which 
15 frequencies represent oblong density variations. 

The energy in such points/areas will directly relate to the number and/or "strength" (grey 
tone" of such density variations. This energy may itself be used to quantify the density 
variations, or the other parts of the power spectrum may be removed and the resulting 
20 power spectrum inversely Fourier transformed to obtain image data representing the ob- 
long density variations. Subsequent to that, these density variations may be quantified as 
mentioned above in relation to the DOG. 

In this respect it should be noted that investigations point to the possibility that certain 
25 types of bone diseases and conditions increase not the number of oblong channels in the 
bone but rather the length or dimensions thereof. Thus, by adapting the dimensions 
sought after using this dimension specific methodology, such variations may be both 
identified and quantified. 

30 Example 3: the DXR BMD technology 

In the following, a description is given of the product Pronosco X-posure System™ that 
incorporates the DXR BMD technique and with reference to the drawings wherein: 

35 Rotate the bone so that it is vertical in the image. 
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An outer edge is found as the pixel path from the top row to the bottom row with highest 
average gradient or highest average curvature in the image (in the relevant side of the 
bone). 

5 When the outer edge is found, the corresponding inner edge is found in at two step pro- 
cedure: 

Firstly, the image rows are shifted horizontally, such that the outer edge pixel in each row 
is aligned with (i.e. right above or right below) the edge pixels in the two neighbouring 
10 rows. The result hereof is called an aligned image. 

Secondly, the aligned image is smoothed in order to reduce the disturbing effect of poten- 
tial resorption spaces/cavities. The smoothing is mainly performed in the vertical direction 
(along the outer edge) using e.g. a standard (flat) local mean mask width dimensions 
15 5mm times 1 mm in the vertical and horizontal directions. Gaussian masks may also be 
applied, and the dimensions can be varied, using e.g. a mask of size 2.5mm times 0.5 
mm or 10mm times 0.5 mm, depending on the image noise level and the dimensions of 
the given bone. 

Thirdly, the inner edge is found in the smoothed aligned image as the pixel path (from the 
20 top row to the bottom row) with highest average grey level value. 

Figure 16 illustrates the cortical thickness t and the bone width W. Shown is a cross sec- 
tion of a cortical bone and the thickness profile of the bone when projected onto the film. 
The projection shows a characteristic shape with local maxima at the boundaries of the 
25 medullar region (at the endosteal surface). 

Figure 17 illustrates the regions used in a Hologic QDR 2000 DEXA apparatus as made 
by Hologic, USA. The MID region starts 15 mm below the lower edge of the radius end 
plate. The bottom of the MID region is defined as 10 mm above the location of the middle 
30 of the 1/3 region. The 1/3 region in turn is defined using information about the length of 
the forearm. See the Hologic Manual in Appendix F for a complete description of the 
Hologic definition of the MID. 

Figure 18 illustrates the regions of interest used in the Pronosco X-posure System™ em- 
35 bodying the technology. Comparing with figure 2, it is seen that the ROIs in the radius and 
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ulna are similar - but not identical - to the regions used in the DEXA apparatus. In the 
metacarpals, the cortical thicknesses are determined on both sides (radial and ulnar 
sides), while in the radius and ulna, only the radial and ulnar sides, respectively, are de- 
termined. 

5 

Theoretical considerations concerning BMP from radioqrammetrv 

To define BMD in the context of e.g. a DEXA machine, one projects the bone onto a 
plane. In this projection, one defines a certain region or area of the bone. BMD is defined 
10 as the mineral mass of the bone projected onto this area, divided by the area itself: 

BMD = Bone Mineral Mass / Area 

The mass used is the ash weight, i.e. the mineral content, which is mainly constituted by 
15 Calcium. The unit of BMD is g/cm 2 or mg/cm 2 . 

When defining the area, care must be taken to understand the projection geometry. In the 
case of X-ray images and fan-beam mode of DEXA machines, there is a certain magnifi- 
cation M of the distances in the image of the bone. This magnification is dictated by the 
20 distance D^from the X-ray emitting spot to the film and the distance D of from the object 
to the film according to the following relation: 

This factor, M, is incorporated in the const-values used for the determination of DXR- 
25 BMD. 

The X-ray images for the X-posure System™ are acquired with a nominal spot-film dis- 
tance D sJ =100 cm. The hand rests directly on the film cassette, and this gives a distance 

D of in the range 1-5 cm. A change of spot-film distance of 10 cm gives a change in M of 

30 0.3% (assuming D of = 3 cm). A change of object-film distance of 2 cm (e.g. by lifting the 

hand from the film) gives a change in M of 2%. Hence it is much more important to control 
the object-film distance than the spot-film distance. This is reflected in the user manual for 
the system. 
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DEXA bone densitometers utilise absorptiometry to estimate the BMD. These devices 
measure the reduction of the intensity of X-rays as they pass though the tissue. This is 
also true for Radiographic Absorptiometry (RA). 

5 In contrast, the Pronosco X-posure System™ does not use the intensities of the image in 
a quantitative manner. This would require the presence of e.g. an aluminum wedge in 
every radiograph, as used in RA-systems. 

Instead, the Pronosco X-posure System™ bases its measurement on radiogrammetry. 
1 0 Radiogrammetry is the measurement of distances on the film. 

The following sections describe the radiogrammetry and the texture analysis. 

Radiogrammetry and bone volume per area (VPA) 

15 

The Pronosco X-posure System™ determines the cortical bone volume using radiogram- 
metry. 

The bone volume per area is denoted VPA. For a cylindrical bone, the following exact 
formula can be derived: 

20 

VPA = * *t*(1-t/W) 

Here Wis the diameter of the bone and t is the cortical thickness, as illustrated in figure 
16. If the bone is not cylindrical, the factor n is replaced by a geometrical factor t 

25 

VPA = f*f *(1-t/W) 

The factor f depends on the shape of the bone. It can be assumed that the cross sectional 
shape of bones for the population covered by the BMD formula (Caucasian women) is in- 
30 variant to a good approximation, i.e. that only the cortical thicknesses and the overall 
sizes vary, while the general shape is constant. With this assumption, the factor f is a 
constant. This assumption is eventually verified by the high correlation between DXR- 
BMD and DEXA-BMD demonstrated in the clinical trials. 

35 The bone volume is used to compute the BMD estimate 
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BMD = c*VPA 

This relation states that the bone volume is proportional to the bone mass. This is based 
on the assumption that cortical bone tissue has a constant mineral content per unit vol- 
5 ume. This assumption is eventually verified in the high correlation between DXR-BMD and 
DEXA-BMD demonstrated in the clinical trials. 

If, however, porosity of the bone is desired taken into account, c could be replaced by c x 
(1-p) where p is a measure of porosity determined from the image data. Otherwise, age 
10 may be introduced, as age and porosity in the dense bone has been found to correlate. 

The measurement region for VPA 

The DEXA apparatus uses the region MID shown in figure 17. 

15 

Both the DXR-BMD and the DEXA-BMD are estimates of the true BMD in the MID region. 
DXR-BMD arrives at this estimate using information from several bone regions. These 
regions are shown in figure 18. 

20 The VPA used in DXR-BMD is computed as an average of VPAs computed in the five 
bones according to the formulae 

VPA = y 2 VPA^uma + V, VPA^t 

25 VPA.net = 1/3 (VPA met2 + VPA ^ + VPA^) 

VPAradiusAilna = % (VPAradius + VPA^na) 

VPA = rc*f*(1-f/M0 

30 

VPA is the bone volume per area and is computed from the average cortical thickness t 
and the outer diameter W of the bone. VPA has the dimension mm. For the metacarpals t 
is the average of the thicknesses in the two sides, while for radius and ulna only the outer 
sides are used. The axial aspects of radius and ulna protrude wing-shaped into the inter- 
35 osseous space to anchor the interosseous membrane. 
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In comparison with traditional radiogrammetry, the Pronosco X-posure System™ utilises 
the automation and computing power of modern PCs to compound measurements from a 
very large number of points along the bones to yield an average cortical thickness for a 
given region of a given bone. The Pronosco X-posure System™ measures the cortical 
5 thickness at 1 1 8 points per centimetre along the axis of a bone, which adds to the repro- 
ducibility and accuracy of the BMD estimate. 

The use of radiogrammetry in the metacarpals to estimate the BMD in the forearm is 
based on the high correlation between BMD in these bones. The same applies to the use 
10 of a slightly different ROI in the radius and ulna as compared to the ROI in the predicate 
device. It should be noted that the ROI's in ulna and radius are positioned identically with 
respect to distances from the 4.0mm landmark. 

In fact, presently, a correlation of 0.90 has been found between the DXR-BMD and DEXA- 
1 5 BMD when using the formula according to the second aspect and as described with re- 
spect to the figures. 

When using the more coarse equation using merely t from the image data, a correlation of 
0.88 is seen between DXR-BMD and DEXA-BMD. 

20 

Example 4: investigation of the effect of certain substances 

Today bone mineral density (BMD) measurements using Dual-energy X-ray Absorpti- 
ometry (DXA) are used to assess the efficacy of anti-osteoporotic treatment with agents 

25 such as Bisphosphonates, calcitonin and HRT. It has been suggested that the changes 
observed in BMD during treatment with these agents only explain part of the observed 
reduction in fracture incidence, and that their effect also might be ascribed to the action on 
the bone micro architecture. A newly developed method termed Digital X-ray Radiogram- 
metry (DXR) may provide more detailed information about treatment effects than obtained 

30 with conventional DXA densitometry. 
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Objectives 

To evaluate the longitudinal changes of different bone parameters measured using DXR 
and DXA technology in selected groups of postmenopausal women, either untreated or in 
5 treatment with hormone replacement therapy (HRT) or Bisphosphonates. 

Methodology 

The study was designed as open comparative, non-randomised observational study. It 
10 was planned that 120 women who had previously participated in a normative reference 
study with a prototype of the Pronosco X-posure System should have their bone status re- 
evaluated using both DXA and DXR. According to treatment status, in the period between 
the initial and the follow-up visit, the women where retrospectively allocated into one of 
the following four groups: an untreated group, a group who received HRT in more than 
15 90% of the study period, a group who received HRT in less than 90% of the study period 
and a group who was treated with Bisphosphonates. 

The following measurements were performed with DXA at the initial visit and at the follow- 
up visit: spine BMD (L2-L4), hip BMD (femoral neck) and forearm BMD (distal). For DXR 
20 the following measurements were performed at the initial visit and at the follow-up visit: 
DXR BMD, cortical thickness of the 2 nd metacarpal (M.T. (2)), porosity and striation. DXA 
BMD was measured using a Norland XR-26 densitometer. DXR BMD, cortical thickness 
of the 2 nd metacarpal, porosity and striation were measured using the Pronosco X-posure 
System. Porosity and striation are measured in arbitrary units on a scale from 0 to 9. 

25 

Based on the calculated mean annual changes the study was analysed with respects to 
changes within the groups as well as between the groups. For within group comparisons a 
paired t-test was used and for between group comparisons an unpaired t-test was used. 
For the within group comparisons SD normalised annual changes were calculated as 

30 (AX 1yea r/SD). 

Ethics 

The study was carried out in accordance with the revised Helsinki Declaration (Somerset 
35 West 1996). The local ethical committee approved the extension of the original protocol 
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and a new patient information and informed consent form before the study was initiated. 
All the participants were given written and oral information about the study before a writ- 
ten informed consent was obtained. 

5 Results 

The baseline demographic characteristics and the DXR BMD values for the subjects en- 
rolled and which can be evaluated are given in table 1 illustrating demographic data and 
observational periods. 

10 

One hundred and twenty-six women were enrolled in the study and 113 were included in 
the analyses. Thirteen women were excluded due to poor quality of the radiographs taken 
at baseline. The main reason for discarding the radiographs was that the hand and fore- 
arm was wrongly centered. Due to the low number of women enrolled in the group who 
15 received HRT in less than 90% of the observational period no further analyses are per- 
formed on this group. In the bisphosphonate group two different agents were used for 
treatment. Eight women got Etidronate and the remaining 3 Alendronate. 
The actual annual changes for the untreated group is shown in table 2 illustrating actual 
annual changes for the untreated group (n=69). 

20 

Expressed in percentages the annual changes for the untreated group is shown in figure 
1. 

The SD normalised annual changes for the untreated group is shown in figure 2. 

25 

The actual annual changes for the group who received HRT in more than 90% of the ob- 
servational period is shown in table 3 illustrating the actual annual changes for the 
"HRT>90%" group (n=28). 

30 Expressed in percentages the annual changes for the "HRT > 90%" group is shown in fig- 
ure 3. 

The SD normalised annual changes for the "HRT > 90%" group is shown in figure 4. 
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The actual annual changes for the bisphosphonate group is shown in table 4 illustrating 
actual annual changes for the bisphosphonate group (n=1 1). 

Expressed in percentages the annual changes for bisphosphonate group is shown in fig- 
5 ure5. 

The SD normalised annual changes for the bisphosphonate group is shown in figure 6. 
The actual annual changes for porosity and striation for all groups are shown in figure 7. 

10 

In order to compare the bisphosphonate treated with the untreated women of comparable 
age, all women from these two groups with an age between 70 and 80 years were ana- 
lysed. In this age matched sub-sample a between groups comparison was made. In table 
5 (illustrating the demographic data and observational periods for the subsample of 
1 5 untreated and bisphosphonate treated women) the demographic baseline data for the two 
sub-samples are shown. 

When comparing the untreated with the bisphosphonate treated groups the actual annual 
changes in DXR BMD, cortical thickness of the 2 nd metacarpal and porosity came out to 
20 be statistically significant difference. Neither spine DXA, femoral neck DXA nor distal 
forearm DXA were statistically significant different The annual changes for the two groups 
and the p-values are given in table 6 illustrating actual annual changes for the subsample 
of untreated and bisphosphonate women. 

25 Expressed in percentages the annual changes for the untreated and the bisphosphonate 
treated group are shown in figure 8. 

The actual annual changes for porosity and striation for the untreated and the bisphos- 
phonate treated groups are shown in figure 9. For porosity a statistically significant reduc- 
30 tion was observed (P=0.046). 

The untreated and the HRT treated group were compared with respect to the actual an- 
nual changes. The demographic data for the two groups are shown in Table 1 . When 
comparing the untreated with the HRT treated group, with respect to the actual annual 
35 changes, cortical thickness of the 2 nd metacarpal was the only parameter showing a stati- 
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cally significant difference. The actual annual changes for the two groups (untreated and 
the "HRT>90%" treated group) and the p-values are given in table 7. 

Expressed in percentages the annual changes for the untreated and the HRT treated 
5 group are shown in figure 10. 

The actual annual changes for porosity and striation for the untreated and the HRT 
treated group are shown in figure 11. 

10 Discussion and Conclusion 

In the present study the Pronosco X-posure System has for the first time been used to 
evaluate longitudinal changes in bones. A total of 126 women were followed for a period 
of two to two-and-a-half year. At baseline and at the end of the study measurements were 
1 5 performed with both DXA and DXR in order to compare the two methods of assessing 
bone status. For the evaluation of the effect of the anti-osteoporotic treatment the study 
participants were retrospectively divided into three groups: an untreated group, a HRT 
group and a bisphosphonate group. 

20 As expected a decrease in BMD of approximately 1 to 1.5% per year was observed in the 
untreated group of women when measured with both DXA and DXR. The annual changes 
in M.T. (2), porosity and striation were also within the expected ranges when data from a 
cross-sectional US normative study are used as reference 1 . Looking at the SD normalised 
annual changes all the BMD measures and the M.T. (2) showed a statistically significant 

25 decrease. 

As a result of the intervention the HRT treated group had a somewhat lower decrease in 
BMD compared to the untreated group. When the SD normalised annual changes were 
compared, hip DXA and forearm DXA showed a statistically significant decrease despite 
30 of treatment, which was in contrast to most of the parameters measured with DXR. One 
exception was striation, which showed a statistically significant increase. M.T. (2) was the 
only measure that was statistically significant difference when the untreated and HRT 
treated group was compared indicating an effect of estrogen on the cortical thickness. 
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For the bisphosphonate treated group an increases in BMD was measured at all sites ex- 
pect for forearm DXA. Also for M.T. (2) a small increase was seen. Both porosity and 
striation showed a decrease overtime. Looking at the SD normalised annual changes po- 
rosity was the only measure, which came out statistically significant. Porosity measured 
5 by DXR is believed to express the porosity of the cortical bone. Base on data from a US 
normative reference study it has been shown that a high level of porosity is associated 
with an increased fracture incidence, which emphasise the possible importance of this 
measure. The study results seems to indicate that treatment with Bisphosphonates have a 
pronounced influence on the cortical porosity. When a direct compassion was made be- 
10 tween the untreated and the bisphosphonate treated group only DXR BMD, M.T. (2) and 
porosity came out statically significant different. 

Results from several clinical studies where pharmacological interventions have been used 
seems to confirm that a simple measurement of bone mass using the DXA technology 

1 5 only provide part of the picture related to the changes which occurs in the bones following 
treatment of osteoporosis. The present clinical study seems to confirm these results by 
indicating that the DXR technology might be more sensitive than DXA to changes in the 
bones caused by anti resorptive treatment with Bisphosphonates and HRT. Further, the 
DXR technology might also be able to distinguish between the individual anti resorptive 

20 agents with respect to their pharmacological effects on the bones. 

In conclusion, the DXR technology in comparison to DXA seems to provide a more com- 
plete picture of the processes in the bones following treatment of osteoporosis, but the 
result form the current observational study must be repeated in a lager population using 
randomised study design. 

25 

Example 5: A system for performing a preferred embodiment of the invention 

In figure 19, a computer system comprising a database 1 is illustrated. The database 1 is 
connected to a means 2 which is used for searching in the database 1 and for outputting 
30 the results of the searches and for receiving searching criteria for use in the search. 

The database 1 will comprise records each identifying a substance having an effect on 
bone and also information relating to one or more parameters affected by the substance. 
The records will also comprise information relating to doses and time intervals, which give 
35 the said effect on the bone. 
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Also, these parameters will be numerical parameters, whereby the information relating to 
a parameter will be a number relating to how much the individual parameters are affected 
by the doses and time intervals. 

5 One substance may be identified by a number of records - where each record relates to 
different doses and/or time intervals between doses. 

The means 2 comprise means for inputting parameters, which should be affected on a 
given bone, which has been determined to be in need of a substance. These values for 
10 these parameters are input, where after the search means identifies substances, which 
affect these parameters. 

It may very well be that a substance will not be present which affects all the given pa- 
rameters in the correct direction - or which affects all parameters to the desired degree. A 

15 substance affecting a number of the parameters may have an adverse effect on one of 
the given parameters - or have an effect on parameters, which are not in the given pa- 
rameters - such as a parameter, which in the bone is normal. In this situation, the search 
means 2 will search for a combination of substances, the combined effect of which has 
the desired overall effect on the given - and preferably also other - parameters of the 

20 bone. 

The search means 2 will finally output the information relating to the substance or sub- 
stances identified. This information will relate to the identity thereof as well as information 
relating to doses and time intervals between doses. 

25 

Naturally, the individual elements of this system - such as the hardware and software, 
may be elements widely in use today, as the overall effect of the system is the identifica- 
tion of suitable parameters and the combination of several substances, if a single sub- 
stance is not identifiable which provides the desired effect. Also, the means 1 and two 
30 may be separated, and a number of search means 2 may be physically distributed and 
connected to the same database 1. In fact, one and the same database 1 may provide the 
information required for search means 2 positioned all over the world. 
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LJ ntreated 


HU']>90% 


HkT<90% 


UisDhos. 


Enrolled subjects (n) 


75 


32 


5 


14 


Evaluable subjects (n) 


69 


28 


5 


11 


Observational period (days) 


868 


842 


904 


804 


Age (years)* 


65.6 (7.8) 


63.6 (7.8) 


56.6 (7.1) 


75.9 (3.3) 


DAK (g/cm2)* 


0.498 (0.06) 


0.548 (0.06) 


0.560 (0.05) 


0.411 (0.06) 


Height (cm)* 


161.8 (5.9) 


163.7 (5.4) 


169.1 (8.0) 


155.8 (5.5) 


Weight (kg)* 


70.8(15.9) 


67.9(10.8) 


71.1 (8.7) 


61.5 (7.1) 



Table 1. 

* Mean (SD) 



5 





Mean 


SD 


DXA Spine (g/cm2) 


-0.0078 


0.0249 


DXA Hip (g/cm2) 


-0.0126 


0.0199 


DXA Forearm (g/cm2) 


-0.0060 


0.0098 


DXR (g/cm2) 


-0.0044 


0.0068 


M.'i'. (2) (mm) 


-0.0199 


0.0390 


Porosity 


-0.0096 


0.4111 


Striation 


0.1113 


0.4859 



Table 2. 



10 Table 3. 





Mean 


Sl> 


DXA Spine (g/cm2) 


-0.0048 


0.0185 


DXA Hip (g/cm2) 


-0.0085 


0.0137 


DXA Forearm (g/cm2) 


-0.0036 


0.0073 


Dxk(g/cm2) 


-0.0016 


0.0072 


M.T. (2) (mm) 


-0.0006 


0.0236 


Forosity 


0.0145 


0.5191 


Striation 


0.1512 


0.3640 



SUBSTITUTE SHEET (RULE 26) 
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Table 4. 





Mean 




DXA Spine (g/cm2) 


0.0154 


0.0294 


DXAHip(g/cm2) 


0.0018 


0.0231 


DXA Forearm (g/cm2) 


-0.0009 


0.0080 


DXk(g/cm2) 


0.0034 


0.0076 


M. 1.(2) (mm) 


0.0024 


0.0113 


Porosity 


-0.3094 


0.3735 


Striation 


-0.0903 


0.4494 



5 Table 5. 





Untreated 


Bisphosphonates 


N 


21 


10 


Observational time (days)* 


836 (125) 


800 (144) 


Age (years)* 


74.3 (2.7) 


76.4(1.6) 


DXk (g/cm2)* 


0.453 (0.056) 


0.411 (0.064) 


Height (cm)* 


162.7 (6.6) 


155.5 (5.8) 


Weight (kg)* 


69.7 (13.3) 


61.0 (7.3) 


Menopausal age (years)* 


48.9 (4.3) 


48.8 (3.0) 



Table 6. 





Untreated 


tfisphosphonates 


P-Values 


N 


21 


10 




DXA Spine (g/cm2)* 


0.0008 (0.0262) 


0.0126 (0.0300) 


0.287 


DXA Hip (g/cm2)* 


-0.0142(0.0211) 


0.0027 (0.0242) 


0.063 


DXA Forearm (g/cm2)* 


-0.0076 (0.0092) 


-0.0011 (0.0084) 


0.07 


DXR (gycm2)* 


-0.0037 (0.0067) 


0.0038 (0.0079) 


0.01 


M. (2) (mm)* 


-0.0154 (0.0320) 


0.0011 (0.0110) 


0.043 


JPorosity* 


-0.0156 (0.3680) 


-0.3169 (0.3928) 


0.046 


Striation* 


0.1660 (0.5093) 


-0.1387 (0.4423) 


0.12 



10 *Mean(SD) 
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Table 7. 





Untreated 


HRT>90% 


IP-Values 


N 


69 


28 




DXA Spine (g/cm2)* 


-0.0078 (0.0249) 


-0.0048 (0.0185) 


0.566 


DXA Hip (g/cm2)* 


-0.0126(0.0199) 


-0.0085 (0.0137) 


0.254 


DXA Forearm (g/cm2)* 


-0.0060(0.001)7) 


-0.0036(0.0073) 


0.253 


DXR(g/cm2)* 


-0.0044(0.0068) 


-0.0016(0.0072) 


0.077 


M.T. (2) (mm)* 


-0.0199(0.0390) 


-0.0006(0.0236) 


0.004 


Porosity* 


0.0096(0.4111) 


0.0145(0.5191) 


0.81 


Striation* 


0.1113(0.4859) 


0.1512(0.3640) 


0.696 
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CLAIMS 

1. A method of determining which substance, of a number of substances having an effect 
on bone, to administer to a vertebrate in need thereof, the method comprising: 

5 - providing image data relating to one or more bones of the vertebrate, 

deriving, from the image data, one or more parameters relating to the bone, 
providing first information relating to the one or more parameters of the one or 

more bones of vertebrates classified as not suffering from bone diseases, 

comparing the first information and the derived one or more parameters and 
10 determining differences there between, 

providing second information relating to each of the number of substances and 

information relating to which of the one or more parameters of bones are affected thereby, 
determining which substance/substances to administer to the vertebrate on the 

basis of the result of the comparison and the second information. 

15 

2. A method according to claim 1, wherein the deriving step comprises deriving, from the 
image data, one or more parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

20 - a variation of a density of cortical and/or trabecular bone of the bone. 

3. A method according to claim 1 or 2, the method being performed on a vertebrate or 
vertebrates classified as suffering from Bone disorders or conditions, such as type I and II 
osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 

25 osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 
diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
such as breast cancer or prostate cancer. 

30 4. A method according to any of claims 1-3, wherein the determination step comprises 
deriving information on suitable doses or intervals of administration. 

5. A method according to claim 1 , wherein the second information further comprises, for 
each substance and for at least one numerical parameter, information relating to whether 
35 the parameter is increased or decreased by the effect of the substance. 
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6. A method according to claim 5, wherein the given parameter(s) comprises the at least 
one numerical parameter and wherein the comparing step also comprises deriving 
information relating to whether a desired effect of the substance is to reduce or increase 
the numerical parameter. 

5 

7. A method according to claim 6, wherein the determining step comprises determining a 
combination of substances providing the desired effect on the bone, where the 
determining step takes into account the information relating to the numerical parameter for 
each of the substances in the combination of substances. 

10 

8. A method according to any of claims 1-7, wherein the determining step comprises: 

providing, in relation to each of the one or more parameters, information relating to 
an uncertainty, such as a statistical uncertainty and/or an uncertainty deriving from the 
providing of the image data or the parameters, in the deriving of the parameter(s) and 
15 - determining, on the basis of that uncertainty, the substance/substances only on 
the basis of parameters) significant in relation to their related uncertainties. 

9. A method according to any of claims 1-8, wherein the determining step comprises 
determining a combination of substances, the effect of which is expected to be altering the 

20 one or more parameters toward those of the first information. 

10. A method according to any of the preceding claims, wherein the step of providing 
second information relating to one or more substances affecting bones and information 
relating to which of the one or more parameters of bones are affected thereby comprises: 

25 - providing first image data relating to at least part of a bone of a vertebrate at a first 
point in time, 

administering one of the one or more substances to the vertebrate after the first 
point in time, 

providing, at a second point in time subsequent to the administering of the 
30 substance, second image data relating to the at least part of the bone, 

deriving, from the first image data, one or more parameters relating to the bone at 
the first point in time, 

deriving, from the second image data, the one or more parameters relating to the 
bone at the second point in time, 
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comparing the one or more parameters relating to the first and second points in 
time, and 

providing, from the result of the comparison, information relating to which of the 
one or more parameters is affected by the substance. 

5 

1 1 . A method according to claim 1 , wherein the deriving steps comprise deriving, from the 
image data, one or more parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

10 - a variation of a density of cortical and/or trabecular bone of the bone. 

12. A method according to claim 10 or 1 1, comprising administering the substance to the 
vertebrate a number of times between the first and second points in time and at 
predetermined doses and time intervals. 

15 

13. A method according to any of claims 10-12, the method further comprising: 

providing third image data of the at least part of the bone and at at least one third 
point in time between the first and second points in time, 

deriving, from the third image data, the one or more parameters, and 
20 - comparing, in the comparing step, the one or more parameters relating to the first, 
second, and third points in time. 

14. A method according to any of claims 1 0-13, wherein: 

the image data are provided relating to at least parts of bones of a number of 
25 vertebrates, 

the substance is provided to each of the vertebrates, 

the one or more parameters are derived relating to the bones of from each of the 
vertebrates, and 

the parameters of the vertebrates are compared in the comparing step. 

30 

15. A method according to any of claims 10-14, wherein the vertebrate or vertebrates 
is/are classified as suffering from Bone disorders or conditions, such as type I and II 
osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 
osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 

35 diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
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diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
such as breast cancer or prostate cancer. 

16. A method according to any of the preceding claims, wherein the step of providing 

5 second information relating to one or more substances affecting bones and information 
relating to which of the one or more parameters of bones are affected thereby comprises: 

providing a first group of vertebrates having received one of the one or more 
substances, 

providing a second group of vertebrates not having received the substance to any 
10 significant degree, 

providing first image data relating to at least part of the bone of each of the first 
group of vertebrates, 

providing second image data relating to the at least part of the bone of each of the 
second group of vertebrates, 
1 5 - deriving, from the first image data, one or more parameters relating to the bones of 
the first group of vertebrates, 

deriving, from the second image data, the one or more parameters relating to the 
bones of the second group of vertebrates, 

comparing the one or more parameters relating to the bones of the first and 
20 second groups of vertebrates and 

providing, from the result of the comparison, information relating to which of the 
one or more parameters are affected by the substance, 

17. A method according to claim 16, wherein the deriving steps comprise deriving, from 
25 the image data, one or more parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

30 18. A method according to claim 16 or 17, wherein the vertebrates of the first group of 
vertebrates have been classified to have or have been brought to have a given disease or 
condition or a disease or condition within a given group of diseases or conditions, and 
wherein the vertebrates of the second group of vertebrates have been classified to not 
have the given disease or condition or a disease or condition within the given group of 

35 diseases or conditions. 
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19. A method according to claim 16 or 17, wherein the vertebrates of the first and second 
groups of vertebrates have been classified to have a given disease or condition or a 
disease or condition within a given group of diseases or conditions. 

5 

20. A method according to claim 19, further comprising: 

providing a third group of vertebrates, wherein the vertebrates of the third group of 
vertebrates have been classified to not have the given disease or condition or a disease 
or condition within the given group of diseases or conditions, 
10 - providing third image data relating to the bone of each of the third group of 
vertebrates, 

deriving, from the third image data, the one or more parameters, and (for hvert 
vertebrat) 

comparing, in the comparing step, the one or more parameters of the first, second, 
15 and third groups of vertebrates. 

21 . A method according to claim 18 or 19 wherein the disease or condition or group of 
diseases or conditions is/are chosen from: Bone disorders or conditions, such as type I 
and II osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 

20 osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 
diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
such as breast cancer of prostate cancer. 

25 22. A method according to any of claims 16-21 , comprising providing the first, second and 
optionally third image data at at least two different points in time between which the first 
group has received at least part of the substance received, and wherein the comparing 
step comprises comparing the parameters derived from the image data obtained at the at 
least two points in time. 

30 

23. A method according to claim any of the preceding claims, wherein the step of 
providing second information relating to one or more substances affecting bones and 
information relating to which of the one or more parameters of bones are affected thereby 
comprises: 
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providing first image data relating to at least part of a bone of a vertebrate at a first 
point in time, 

administering one of the one or more substances to the vertebrate after the first 
point in time, 

5 - providing, at a second point in time subsequent to the providing of the substance, 
second image data relating to the at least part of the bone, 

deriving, from the first image data, one or more parameters relating to the bone at 
the first point in time, 

deriving, from the second image data, the one or more parameters relating to the 
10 bone at the second point in time, 

comparing the one or more parameters relating to the first and second points in 
time, and 

determining, from the result of the comparison, whether the substance has any 
effect on the bone, and 
15 - providing, if the substance has an effect on the bone, information that the 

substance has an effect on the bone and which of the one or more parameters is affected 
thereby. 

24. A method according to claim 23, wherein the deriving steps comprise deriving, from 
20 the image data, one or more parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

25 25. A method according to claim 23 or 24, further comprising administering the substance 
to the vertebrate a number of times between the first and second points in time and at 
predetermined doses and time intervals. 

26. A method according to claim 23, 24, or 25, the method further comprising: 
30 - providing third image data of the at least part of the bone and at at least one third 

point in time between the first and second points in time, 

deriving, from the third image data, the one or more parameters, and 
comparing, in the comparing step, the one or more parameters relating to the first, 

second, and third points in time. 

35 
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27. A method according to any of claims 23-26, wherein: 

the image data are provided of at least parts of bones of a number of vertebrates, 

the substance is provided to each of the number of vertebrates, 

the one or more parameters are derived from each of the number of vertebrates, 

5 and 

the parameters of the vertebrates are compared in the comparing step. 

28. A method according to any of the preceding claims, wherein the vertebrate or 
vertebrates is/are classified as suffering from Bone disorders or conditions, such as type I 

10 and II osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 
osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 
diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
such as breast cancer or prostate cancer. 

15 

29. A method according to any of the claims, wherein the step of providing second 
information relating to one or more substances affecting bones and information relating to 
which of the one or more parameters of bones are affected thereby comprises: 

providing a first group of vertebrates having received one of the one or more 
20 substances, 

providing a second group of vertebrates not having received the substance to any 
significant degree, 

providing first image data relating to at least part of the bone of each vertebrate of 
the first group of vertebrates, 
25 - providing second image data relating to the at least part of the bone of each 
vertebrate of the second group of vertebrates, 

deriving, from the first image data, one or more parameters relating to the bones of 
the first group of vertebrates, 

deriving, from the second image data, the one or more parameters relating to the 
30 bones of the second group of vertebrates, 

comparing the one or more parameters relating to the first and second groups of 
vertebrates, 

determining, from the result of the comparison, whether the substance has any 
effect on the bones of the vertebrates of the first group, and 
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providing, if the substance has an effect on the bone, information that the 
substance has an effect on the bone and which of the one or more parameters is affected 
thereby. 

5 30. A method according to claim 29, wherein the deriving steps comprise deriving, from 
the image data, one or more parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

10 

31 . A method according to claim 29 or 30, wherein the vertebrates of the first group of 
vertebrates have been classified to have a given disease or condition or a disease or 
condition within a given group of diseases or conditions, and wherein the vertebrates of 
the second group of vertebrates have been classified to not have the given disease or 

15 condition or a disease or condition within the given group of diseases or conditions. 

32. A method according to claim 29, wherein the vertebrates of the first and second 
groups of vertebrates have been classified to have a given disease or condition or a 
disease or condition within a given group of diseases or conditions. 

20 

33. A method according to claim 32, further comprising: 

providing a third group of vertebrates, wherein the vertebrates of the third group of 
vertebrates have been classified to not have the given disease or condition or a disease 
or condition within the given group of diseases or conditions, 
25 - providing third image data relating to the bone of each of the third group of 
vertebrates, 

deriving, from the third image data, the one or more parameters, and 
comparing, in the comparing step, the one or more parameters of the first, second, 
and third groups of vertebrates. 

30 

34. A method according to claim 31 or 32, wherein the disease or condition or group of 
diseases or conditions is/are chosen from Bone disorders or conditions, such as type I 
and II osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 
osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 

35 diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
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diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
such as breast cancer or prostate cancer. 

35. A method according to any of claims 29-34, further comprising providing the first, 
5 second and optionally third sets of image data at at least two different points in time 
between which the first group has received at least part of the substance received, and 
wherein the comparing step comprises comparing the image data obtained at the at least 
two points in time. 

10 36. A method of differentiating between different types of bone lesions/defects/conditions, 
the method comprising: 

a) providing image data on a bone in at least one vertebrate, which is suspected to suffer 
from a first bone disease or condition, 

b) providing image data on a bone in at least one vertebrate, which is suspected to suffer 
15 from a second bone disease or condition 

c) deriving one or more parameters relating to the bones based on said image data 
provided in step a) and b), 

d) providing information relating to said one or more parameters of bones of vertebrates 
not suffering from the bone disease or condition, 

20 e) comparing the one or more parameters obtained in step b) with the one or more 
parameters obtained in step c) to determine differences in said parameters, and 
f) identifying, on the basis of the information provided in step d) and the result of the 
comparison of step e), one or more parameters relating to said bones of said vertebrates 
and being correlated to one or each of the first and second bone diseases or conditions. 

25 

37. A method according to claim 36, wherein the deriving step comprises deriving, from 
the image data, one or more parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

30 - a variation of a density of cortical and/or trabecular bone of the bone. 

38. A method of determining the effect of a substance on a bone of a vertebrate, the 
method comprising: 

providing first image data relating to at least part of the bone at a first point in time, 
35 - administering the substance to the vertebrate after the first point in time, 
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providing, at a second point in time subsequent to the administering of the 
substance, second image data relating to the at least part of the bone, 

deriving, from the first image data, one or more parameters relating to the bone at 
the first point in time, 

5 - deriving, from the second image data, the one or more parameters relating to the 
bone at the second point in time, 

comparing the one or more parameters relating to the first and second points in 
time, and 

determining, from the result of the comparison, the effect of the substance on the 

10 bone, 

wherein the deriving steps comprise deriving, from the image data, one or more 
parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

15 and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

39. A method according to claim 38, comprising administering the substance to the 
vertebrate a number of times between the first and second points in time and at 

20 predetermined doses and time intervals. 

40. A method according to claim 38 or 39, the method further comprising: 

providing third image data of the at least part of the bone and at at least one third 
point in time between the first and second points in time, 
25 - deriving, from the third image data, the one or more parameters, and 

comparing, in the comparing step, the one or more parameters relating to the first, 
second, and third points in time. 

41 . A method according to any of claims 38-40, wherein: 

30 - the image data are provided relating to at least parts of bones of a number of 
vertebrates, 

the substance is provided to each of the vertebrates, 

the one or more parameters are derived from each of the vertebrates, and 

the parameters of the vertebrates are compared in the comparing step 
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in order to provide information relating to an overall effect of the substance on bone 
statistical analysis of the parameters of the vertebrates. 

42. A method according to any of claims 38-41 , wherein the vertebrate or vertebrates 
5 is/are classified as suffering from Bone disorders or conditions, such as type I and II 

osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 
osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 
diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
10 such as breast cancer or prostate cancer. 

43. A method of determining the effect of a substance on a bone of a vertebrate, the 
method comprising: 

providing a first group of vertebrates having received the substance, 
15 - providing a second group of vertebrates not having received the substance to any 
significant degree, 

providing first image data relating to at least part of the bone of each of the first 
group of vertebrates, 

providing second image data relating to the at least part of the bone of each of the 
20 second group of vertebrates, 

deriving, from the first image data, one or more parameters relating to the bones of 
the first group of vertebrates, 

deriving, from the second image data, the one or more parameters relating to the 
bones of the second group of vertebrates, 
25 - comparing the one or more parameters relating to the bones of the first and 
second groups of vertebrates, and 

determining, from the result of the comparison, the effect of the substance on the 
bones of the vertebrates of the first group, 

30 wherein the deriving steps comprise deriving, from the image data, one or more 
parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

35 



WO 00/72216 



65 



PCT/DKOO/00273 



44. A method according to claim 43, wherein the vertebrates of the first group of 
vertebrates have been classified to have or have been brought to have a given disease or 
condition or a disease or condition within a given group of diseases or conditions, and 
wherein the vertebrates of the second group of vertebrates have been classified to not 

5 have the given disease or condition or a disease or condition within the given group of 
diseases or conditions. 

45. A method according to claim 43, wherein the vertebrates of the first and second 
groups of vertebrates have been classified to have a given disease or condition or a 

10 disease or condition within a given group of diseases or conditions. 

46. A method according to claim 45, further comprising: 

providing a third group of vertebrates, and wherein the vertebrates of the third 
group of vertebrates have been classified to not have the given disease or condition or a 
1 5 disease or condition within the given group of diseases or conditions, 

providing third image data relating to the bone of each of the third group of 
vertebrates, 

deriving, from the third image data, the one or more parameters, and 
comparing, in the comparing step, the one or more parameters of the first, second, 
20 and third groups of vertebrates. 

47. A method according to claim 44 or 45 wherein the disease or condition or group of 
diseases or conditions is/are chosen from: Bone disorders or conditions, such as type I 
and II osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 

25 osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 
diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
such as breast cancer of prostate cancer. 

30 48. A method according to any of claims 43-47, comprising providing the first, second and 
optionally third image data at at least two different points in time between which the first 
group has received at least part of the substance received, and wherein the comparing 
step comprises comparing the parameters derived from the image data obtained at the at 
least two points in time. 



35 
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49. A method of screening a test substance to identify whether the substance has an 
effect on a bone of a vertebrate, the method comprising: 

providing first image data relating to at least part of the bone at a first point in time, 

administering the substance to the vertebrate after the first point in time, 
5 - providing, at a second point in time subsequent to the providing of the drug, 
second image data relating to the at least part of the bone, 

deriving, from the first image data, one or more parameters relating to the bone at 
the first point in time, 

deriving, from the second image data, the one or more parameters relating to the 
1 0 bone at the second point in time, 

comparing the one or more parameters relating to the first and second points in 
time, and 

determining, from the result of the comparison, whether the substance has any 
effect on the bone, 

15 

wherein the deriving steps comprise deriving, from the image data, one or more 
parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

20 - a variation of a density of cortical and/or trabecular bone of the bone. 

50. A method according to claim 49, further comprising administering the substance to the 
vertebrate a number of times between the first and second points in time and at 
predetermined doses and time intervals. 

25 

51. A method according to claim 49 or 50, the method further comprising providing third 
image data of the at least part of the bone and at at least one third point in time between 
the first and second points in time, deriving, from the third image data, the one or more 
parameters, and comparing, in the comparing step, the one or more parameters relating 

30 to the first, second, and third points in time. 

52. A method according to any of claims 49-51, wherein the image data are provided of at 
least parts of bones of a number of vertebrates, where the substance is provided to the 
number of vertebrates, where the one or more parameters are derived from each of the 

35 vertebrates, and wherein the parameters of the vertebrates are compared in the 
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comparing step in order to provide information relating to an overall effect of the 
substance on bone. 

53. A method according to any of claims 49-52, wherein the vertebrate or vertebrates 
5 is/are classified as suffering from Bone disorders or conditions, such as type I and II 

osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 
osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 
diabetes (IODM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
diseases such as osteogenesis imperfecta, Paget's bone disease, myeloma or cancer, 
10 such as breast cancer or prostate cancer. 

54. A method of screening a test substance to identify whether the substance has an 
effect on a bone of a vertebrate, the method comprising: 

providing a first group of vertebrates having received the substance, 
15 - providing a second group of vertebrates not having received the substance to any 
significant degree, 

providing first image data relating to at least part of the bone of each vertebrate of 
the first group of vertebrates, 

providing second image data relating to the at least part of the bone of each 
20 vertebrate of the second group of vertebrates, 

deriving, from the first image data, one or more parameters relating to the bones of 
the first group of vertebrates, 

deriving, from the second image data, the one or more parameters relating to the 
bones of the second group of vertebrates, 
25 - comparing the one or more parameters relating to the first and second groups of 
vertebrates and 

determining, from the result of the comparison, whether the substance has any 
effect on the bones of the vertebrates of the first group 

30 wherein the deriving steps comprise deriving, from the image data, one or more 
parameters relating to: 

one or more physical distances in the bone or the image data relating to the bone 

and/or 

a variation of a density of cortical and/or trabecular bone of the bone. 

35 
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55. A method according to claim 54, wherein the vertebrates of the first group of 
vertebrates have been classified to have a given disease or condition or a disease or 
condition within a given group of diseases or conditions, and wherein the vertebrates of 
the second group of vertebrates have been classified to not have the given disease or 

5 condition or a disease or condition within the given group of diseases or conditions. 

56. A method according to claim 54, wherein the vertebrates of the first and second 
groups of vertebrates have been classified to have a given disease or condition or a 
disease or condition within a given group of diseases or conditions. 

10 

57. A method according to claim 56, further comprising: 

providing a third group of vertebrates, and wherein the vertebrates of the third 
group of vertebrates have been classified to not have the given disease or condition or a 
disease or condition within the given group of diseases or conditions, 
1 5 - providing third image data relating to the bone of each of the third group of 
vertebrates, 

deriving, from the third image data, the one or more parameters, and 
comparing, in the comparing step, the one or more parameters of the first, second, 
and third groups of vertebrates. 

20 

58. A method according to claim 55 or 56 wherein the disease or condition or group of 
diseases or conditions is/are chosen from Bone disorders or conditions, such as type I 
and II osteoporosis, glucocroticoid induced osteoporosis and other types of drug induced 
osteoporosis, osteomalacia, hyper- and hypoparathyroidism, hyper- and hypothyroidism, 

25 diabetes (IDDM), hypopituitarism, osteoarthritis, rheumatoid arthritis, genetic bone 
diseases such as osteogenesis imperfecta, Pagefs bone disease, myeloma or cancer, 
such as breast cancer or prostate cancer. 

59. A method according to any of claims 54-58, further comprising providing the first, 
30 second and optionally third sets of image data at at least two different points in time 

between which the first group has received at least part of the substance received, and 
wherein the comparing step comprises comparing the image data obtained at the at least 
two points in time. 
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60. A method according to any of the preceding claims, wherein each vertebrate is a 
laboratory or test animal, and wherein the method comprises finally sacrificing the animal. 

61 . A method according to any of the preceding claims, wherein the image data are 

5 obtained by a non-invasive method, such as by exposing the bone or bones to radiation, 
such as X-rays, performing a CT scanning, a NMR scanning or exposing the bone or 
bones to ultra sound. 

62. A method according to claim 61, wherein the image data are two-dimensional image 
10 data relating to at least part of the bone or bones. 

63. A method according to any of the preceding claims, wherein the image data comprise 
information relating to the cortical bone of at least a part of the bone, and wherein one of 
the one or more parameters relates to a thickness, t, of the cortical bone. 

15 

64. A method according to any of the preceding claims, wherein the image data comprise 
information relating to the cortical bone of at least a part of the bone, and wherein one of 
the one or more parameters relates to a width, w, of the bone. 

20 65. A method according to claim 63, wherein the bone is a radius, and wherein the t value 
or values is/are determined on a radial side of the radius. 

66. A method according to claim 63, wherein the bone is an ulna, and wherein the t value 
or values is/are determined on an ulnar side of the ulna. 

25 

67. A method according to claims 63 and 64, wherein t and w are determined along a 
single line extending in a direction perpendicular to the bone. 

68. A method according to claim 67, wherein pairs of (t,w) are determined for a plurality of 
30 lines extending in the direction perpendicular to the bone and being positioned at different 

positions along the longitudinal direction of the bone, and wherein at least one of the one 
or more parameters is determined on the basis of pairs of (t,w) corresponding to the 
individual lines. 
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69. A method according to claim 68, wherein the at least one parameter is determined on 
the basis of mean values of the t and w values corresponding to the individual lines. 

70. A method according to claim 68 or 69, wherein t and w values are determined for each 
5 of more than 10, such as more than 20, such as more than 40, such as more than 60, 

such as more than 80, such as more than 100, such as more than 110 lines per cm. of the 
bone within a predetermined longitudinal part of the bone. 

71. A method according to any of the preceding claims, wherein one or more parameters 
10 relate to a variation value and wherein at least one variation value is related to a cortical 

porosity of the bone. 

72. A method according to any of the preceding claims, wherein one or more parameters 
relate to a variation value and the at least one variation value is determined by a method 

15 comprising: 

determining, within the image data, first parts thereof relating to cortical bone and 
second parts positioned within the first parts, the second parts not relating to cortical 
bone. 

20 73. A method according to any of the preceding claims, wherein the image data are 
derived from a digitised, two-dimensional image of the at least part of the bone, and 
wherein one or more parameters relate to a variation value and at least one variation 
value is derived by determining one or more first areas in the image data relating to or 
representing cortical tissue and one or more second areas positioned within the first 

25 area(s), each of the second area(s) representing a local deviation in grey value. 

74. A method according to claim 73, wherein the local deviation in grey value is a grey 
value being lower than the grey values of any surrounding areas of the first area(s), such 
as lower than the grey values of the surrounding areas minus a predetermined value. 

30 

75. A method according to any of the preceding claims, wherein one or more parameters 
relate to a variation value and at least one variation value is determined by, in the image 
data: 

determining an outer cortical edge of the bone and denoting pixels thereon as 
35 cortical pixels, 
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repeating: 

for each current, cortical pixel, denote any of its neighbouring pixels having 
a grey value higher than that of itself as a cortical pixel, 

until no non-cortical pixels in the neighbourhood of any current cortical 
5 pixel have a grey level value that is greater than the grey level value of the given current 
cortical pixel, thereby defining an inner cortical edge as the inner (endosteal) boundary of 
the final set of cortical pixels, 

denoting non-cortical image data pixels within the outer and inner cortical edges 
as porosity pixels, and 

10 - deriving the at least one variation value on the basis of the cortical pixels and the 
porosity pixels. 

76. A method according to claim 75, wherein the deriving step comprises deriving the at 
least one variation value on the basis on a number of porosity pixels compared to a 

15 number of cortical pixels in the image data, 

77. A method according to claim 75, wherein the deriving step comprises deriving the at 
least one variation value on the basis of a number, size and/or shape of groups of porosity 
pixels. 

20 

78. A method according to claim 77, wherein the groups of porosity pixels are isolated 
groups of porosity pixels being surrounded by cortical pixels. 

79. A method according to any of the preceding claims, wherein one or more parameters 
25 relate to a variation value and at least one variation value is determined from oblong 

density variations in the image data, such as of parts of the image data relating to cortical 
tissue. 

80. A method according to claim 79, wherein the at least one value is determined from 
30 oblong density variations extending at least substantially in a longitudinal direction of the 

bone. 

81. A method according to claim 79 or 80, wherein the oblong density variations have a 
length to width ratio of at least 1 VzA , such as at least 2:1 , such as at least 2V£:1 , 
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preferably at least 3:1, such as at least 3.5:1, preferably at least 4:1, such as at least 
4.5:1. 

82. A method according to any of claims 79-81 , wherein the oblong density variations 
5 have a length to width ratio in the interval 1 J4-8: 1 , such as 1 .7-6: 1 , preferably 2-4: 1 . 

83. A method according to any of claims 79-82, wherein the oblong density variations in 
the image data relate to oblong density variations, such as cavities, tunnels, channels, 
valleys, grooves, or the like, in the cortical bone. 

10 

84. A method according to claim 83, wherein the density variations of the bone have a 
length in the interval of 0.4-15 mm, such as in the interval 0.5-12 mm, preferably in the 
interval 0.8-10 mm, such as in the interval 1-8 mm, such as on the order of 2 mm, and a 
width in the interval 0.1-6 mm, such as in the interval 0.13-5 mm, preferably in the interval 

15 0.15-3 mm, such as in the interval 0.17-2 mm, preferably in the interval 0.18-1 mm, such 
as on the order of 0.2-0.4 mm. 

85. A method according to any of claims 79-84, wherein the oblong density variations of 
image data have a length in the interval of 0.4-15 mm, such as in the interval 0.5-12 mm, 

20 preferably in the interval 0.8-10 mm, such as in the interval 1-8 mm, such as on the order 
of 2 mm, and a width in the interval 0.1-6 mm, such as in the interval 0.13-5 mm, prefer- 
ably in the interval 0.15-3 mm, such as in the interval 0.17-2 mm, preferably in the interval 
0.18-1 mm, such as on the order of 0.2-0.4 mm. 

25 86. A method according to any of claims 79-85, wherein the determined variation value is 
determined from a number of oblong density variations in the image data, a mean width 
thereof, and/or a mean length thereof. 

87. A method according to any of claims 79-86, wherein the determined variation value is 
30 determined from a grey value difference between the oblong density variations and other 

parts of the image data. 

88. A method according to any of claims 79-87, wherein the variation value is determined 
by: 

35 - obtaining a power spectrum relating to the image data, 
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identifying parts thereof relating to frequencies or frequency intervals 
corresponding to predetermined dimensions or dimension intervals of the oblong 
variations by: 

a) estimating an energy at or within the frequencies or frequency 
5 intervals, and/or 

b) removing at least substantially all other parts of the power 

spectrum, and 

inversely Fourier transforming the resulting power spectrum, and 
determining a variation of the grey levels of the resulting image. 

10 

89. A method according to any of claims 79-87, wherein the variation value is determined 
by filtering the image data with an oblong Gaussian Kernel and subsequently determining 
a variation of grey levels of the resulting image. 

15 90. A method according to claim 89, wherein the variation value is determined by filtering 
each of the image data and the resulting image with an oblong Gaussian Kernel, the two 
oblong Gaussian kernels having different dimensions, and subsequently subtracting the 
two resulting image data and determining a variation of grey levels of the resulting image. 

20 91. A method according to claim 89 or 90, wherein the oblong kemel(s) has/have a length 
to width ratio of at least 114:1 , such as at least 2:1 , such as at least 214:1 , preferably at 
least 3:1, such as at least 3.5:1, preferably at least 4:1, such as at least 4.5:1. 

92. A method according to any of claims 89-91 , wherein the oblong kernel(s) has/have a 
25 length to width ratio in the interval 114-8:1, such as 1.7-6:1, preferably 2-4:1. 

93. A method according to any of claims 71-92, wherein the or one oblong kernel has a 
length in the interval of 0.4-15 mm, such as in the interval 0.5-12 mm, preferably in the 
interval 0.8-10 mm, such as in the interval 1-8 mm, such as on the order of 2 mm, and a 

30 width in the interval 0.1-6 mm, such as in the interval 0.13-5 mm, preferably in the interval 
0.15-3 mm, such as in the interval 0.17-2 mm, preferably in the interval 0.18-1 mm, such 
as on the order of 0.2-0.4 mm. 

94. A method according to claim 90, wherein one oblong kernel has a length in the 

35 interval of 1-15 mm, such as in the interval 1 .3-10 mm, preferably in the interval 114-6 mm, 



WO 00/72216 



74 



PCT/DK00/00273 



such as in the interval 2-4 mm, such as on the order of 3 mm, and a width in the interval 
0.4-6 mm, such as in the interval V£-5 1 /4 mm, preferably in the interval 0.7-5 mm, such as 
in the interval 1-2% mm, preferably in the interval 1.2-2 mm, such as on the order of 1!4 
mm. 

5 

95. A method according to any of claims 79-94, wherein the at least one variation value is 
determined by: 

providing a pre-determined template of the oblong density variations, 
template matching parts of the image data with the predetermined shape so as to 
10 identify oblong density variations, 

determining the variation value from: 

the number of identified oblong density variations, 
a grey value difference between the identified oblong density variations 
and other parts of the image data, or 
15 - an average length and/or width of the identified oblong density variations. 

96. A method according to any of the preceding claims, wherein the image data 
comprising information relating to the trabecular structure of at least a part of a bone of 
the vertebrate, and wherein one or more parameters relate to a variation value and at 

20 least one variation value is derived by one or more of the following methods: 

obtaining an estimate of the parametric estimate of the power spectrum of the 
image data and extracting at least one of the one or more parameters as one relating to 
the energy distribution of the parametric estimate, 

obtaining, on the basis of image data on which a Fourier method has been used to 
25 emphasise the information in the image data relating to the trabecular structure, an 
estimate of a grey-level co-occurrence matrix and extracting at least one of the one or 
more parameters on the basis of the estimated co-occurrence matrix, 

obtaining an estimate of the projected trabecular pattern of the image data by 
using a Fourier method to emphasise the information in the image data relating to the 
30 trabecular structure and subjecting the manipulated image data to a morphological 

operation, and extracting at least one of the one or more parameters as one relating to the 
trabecular structure from the estimated projected trabecular pattern, and 
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obtaining, on the basis of a frequency analysis of the image data, at least one of 
the one or more parameters as one relating to the periodicity of the trabecular structure of 
the part of the bone. 

5 97. A method according to any of the preceding claims, wherein the image data 

comprising information relating to the trabecular structure of at least a part of a bone of 
the vertebrate, and wherein one or more parameters relate to a variation value and at 
least one variation value is derived by: 

determining values reflecting the projected trabecular density in the image data, 
10 caused by the X-ray attenuating properties of cancellous bone in the part of the bone, for 
each of a number of locations or areas in the image data, and 

deriving the at least one value from the variation of the determined PTD-values, 
preferably in the longitudinal direction of the bone. 

15 98. A method according to any of the preceding claims, further including determining, from 
the result of the comparison, which of the one or more parameters is/are effected by the 
substance. 

99. A method according to claim 98, wherein a statistical uncertainty is derived in relation 
20 to each of the one or more parameters, and wherein the determined parameter(s) is/are 

that/those having an effect exceeding the uncertainty(ies) relating to the parameters) 

100. A method according to any of the preceding claims, wherein the substance is a 
substance known to or expected to have an effect on bones. 

25 

101. A method according to claim 100, wherein the substance is known to or expected to 
be: 

an antiresorptive agent, such as an agent expected to act through a direct or 
indirect inhibition of the osteoclasts, such as: 
30 - Bisphosphonates, such as Etidronate, Clodronate, Pamidronate, 

Alendronate, Risedronate, Ibandronate, Zoledronate, Tiludronate, Incadronate, 
Neridronate, Olpadronate, EB-1053, or YH 529, 

Hormone Replacement Therapy, such as Estrogens, Estrogens with 
progestogens, Tibolone, or Trimegestone, 
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Selective Estrogen Receptor Modulator (SERM), such as Raloxifene, 
Droloxifene, Idoxifene, Levormeloxifene, or Tamoxifen, or 

Calcitonin, such as Salmon calcitonin, Human calcitonin, Calcitonin gene- 
related peptide, 

5 - an agent with Heterogeneous Effect, such as Stanozolol, Oxandrolone, 

Nandrolone, Dihydrotachysterol and other vitamin D derivatives, Alfacalcidiol, Calcifero, 
Calcitriol, Cholecalciferol, ST-630, Calcium, Calcium with vitamin D, Ipriflavone and other 
isoflavones, or Strontium, 

a bone forming agent, such as an agent expected to act through a direct or indirect 
10 stimulation of the osteoblasts through a myogenic mechanism, such as Sodium fluoride, 
Monofluorophosphate, Parathyroid hormone^, Parathyroid hormone^, Growth 
hormone, Growth hormone releasing compounds, IGF-1, IGF-BP3, Osteogenic protein-1, 
or Bone morphogenetic protein-2. 

15 1 02. A method according to any of the preceding claims, wherein the vertebrate is a 
human, a horse, a large ape, a great ape, an anthropoid ape, a pig, a cow, a rat, a rabbit, 
a dog, a cat. 

103. A method according to any of the preceding claims, wherein the bone is taken from 
20 the group consisting of radius, femur, corpus vertebrae (L1 , L2, L3, L4, L5, T1 , T2, T3, T4, 

T5 f T6, T7, T8, T9, T10, T11, T12, C1 t C2, C3, C4, C5, C6, C7), calcaneus, talus, os 
carpi, metatars, metacarpi, phalanges, tibia, fibula, patella, ulna, humerus, mandible, 
clavicula, scapula, os coxae, os naviculare, os cuboideum, os cuneiform I, os cuneiform II, 
os cuneiform III, os sacrum, os coccygis. 

25 

104. A method according to any of the preceding claims, wherein the image data relate to 
information relating to bone tissue through a thickness of the bone. 

105. A system for identifying one or more substances for treatment of bone 
30 diseases/conditions, the system comprising: 

a data storage means for holding identifiers for a number of substances and, for 
each substance, data relating to one or more bone parameters affected by the substance 
when administered to a vertebrate, 
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search means for comparing one or more given bone parameters to parameters in 
the data storage means in order to identify one or more substances affecting the one or 
more given parameters, 

means for outputting information relating to the identified substance(s). 

5 

106. A system according to claim 105, wherein the storage means are adapted to hold 
data relating to one or more bone parameters for each substance, each parameter 
relating to: 

one or more physical distances in the bone or the image data relating to the bone 

10 and/or 

a variation of a density of cortical and/or trabecular bone of the bone, and 
wherein each of the one or more given parameters is one of the plurality of parameters. 

107. A system according to claim 105 or 106, wherein the search means are adapted to 
15 identify a combination of substances affecting the one or more given parameters of the 

bone. 

108. A system according to any of claims 105-107, wherein the data storage is adapted to 
hold data relating to at least one numerical parameter for each substance as well as 

20 information relating to whether the at least one numerical parameter is increased or 
decreased by the effect of the substance. 

109. A system according to claim 108, wherein the given parameters) comprise(s) at 
least one of the at least one numerical parameter, and wherein the search means further 

25 comprise means for receiving information relating to whether a desired effect of the 
substance is to reduce or increase the at least one numerical parameter. 

110. A system according to claim 108 or 109, wherein the search means comprise means 
for determining a combination of substances providing a desired effect on the bone, 

30 where the search means are adapted to take into account the information relating to the 
numerical parameter. 

1 1 1 . A system according to any of claims 105-1 10, wherein the data storage further 
comprises, for each substance, information relating to one or more recommendable doses 

35 to be administered to a person. 
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112. A system according to any of claims 105-111, wherein the data storage further 
comprises, for each substance, information relating to recommendable time periods 
between doses to be administered to a person. 

5 1 13. A system according to any of claims 105-1 12, wherein the data storage further 
comprises, for each substance, information relating to any known side effects of the 
substance when administered to a person. 

1 14. An apparatus for providing information relating to the effect of a substance on a 
10 vertebrate, the apparatus comprising: 

means for providing image data relating to one or more bones of the vertebrate, 
means for deriving, from the image data, one or more parameters relating to the 

bone, 

means for providing first information relating to the one or more parameters of the 
1 5 one or more bones of vertebrates classified as not suffering from bone diseases, 

means for comparing the first information and the derived one or more parameters 
and determining differences there between, 

means for providing second information relating to each of the number of 
substances and information relating to which of the one or more parameters of bones are 
20 affected thereby, 

means for determining which substance/substances to administer to the vertebrate 
on the basis of the result of the comparison and the second information. 

1 15. An apparatus for determining the effect of a substance on a bone of a vertebrate, the 
25 apparatus comprising: 

means for providing first image data relating to at least part of the bone at a first 
point in time, (preferably - prior to) 

means for providing, at a second point in time, second image data relating to the 
at least part of the bone, (preferably - subsequent) 
30 - means for deriving, from the first image data, one or more parameters relating to 
the bone at the first point in time, 

means for deriving, from the second image data, the one or more parameters 
relating to the bone at the second point in time, 
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the deriving means being adapted to derive, from the image data, one or more 
parameters relating to: 

one or more physical distances in the bone or the image data relating to 
the bone and/or 

5 a variation of a density of cortical and/or trabecular bone of the bone 

means for comparing the one or more parameters relating to the first and second 
points in time, and 

means for determining, from the result of the comparison, the effect of the 
substance on the bone. 
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